»V  s  C* 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


X 


CONTENTS 

Solid  State  Session  P-1 

-  Page 

Initial  Susceptibility  Investigation  of  Magnetic  Transitions  in 
Several  Rare  Earth  Metals:  Thermal  Hysteresis  in  Ferro¬ 
magnetic  Transitions 

F.  J.  Jelinek,  E.  D.  Hill  and  B.  C.  Gerstein .  1 

Distribution  of  Magnetic  Moment  Density  in  Tb 

O.  Steinsvoll,  G.  Shirane,  R.  Nathans .  15 

Magnetostriction  in  the  Heavy  Rare  Earth  Metals 

B.  F.  DeSavage  and  A.  E.  Clark . . .  23 

Observation  of  Optical  Absorption  Bands  Associated  with 
Magnetic  Order  in  Rare-Earth  Metals 

C.  Chr.  Schuler . . .  35 

Infrared  Absorption  Structure  in  Rare  Earth  Metals:  Relationship 
to  Spin  Arrangement  and  Band  Structure 

B.  R.  Cooper  and  R.  W.  Redington .  45 

Electron-Spinwave  Scattering  in  Gadolinium 

B.  Luthi  and  H.  Rohrer .  55 

Thermal  Conductivity  and  Electrical  Resistivity  of  Samarium 
at  Low  Temperatures 

Sigurds  Arajs  and  G.  R.  Dunmyre .  63 

F.esults  of  Augmented  Plane  Wave  Calculations  of  the  Band 
Structure  of  Cerium  Metal 

J.  T.  Waber  and  A.  C.  Switendick .  75 

Theoretical  Results  for- Positron  Annihilation  in  Rare  Earth 
Metals 

T.  L.  Loucks .  89 

High  Pressure  Polymorphism  and  Magnetic  Ordering  in  Some 
Sm-Type  Rare  Earth  Alloys  and  Samarium 

A.  Jayaraman,  R.  C.  Sherwood,  and  E.  Corenzwit .  91 

Valence  Bond  Formation  in  the  Close-Packed  Rare  Earth 
Metals 

Forrest  L.  Carter 


103 
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THERMAL  HYSTERESIS  IN  FERROMAGNETIC 
TRANSITIONS* 

F.  J.  JELINEK,  E.  D.  HILL  and  B.  C.  GERSTEIN 

Institute  for  Atomic  Research  and  Department  of  Physical  Chemistry 

Iowa  State  University,  Ames,  Iowa 
(Work  performed  in  the  Ames  Laboratory  of  the  U.  S.  A.  E.  C.  ) 

ABSTRACT 

Thermal  hysteresis  of  magnetic  transitions  in  polycrystalline 
Gd,  Tb,  Dy,  Ho,  Tm  and  Sm  have  been  investigated  by  an  a.  c.  mutual 
inductance  magnetic  susceptibility  technique.  A  thermal  hysteresis 
effect  was  observed  in  the  ferromagnetic  transitions  of  Tb,  Dy,  Ho  and 
Tm.  No  thermal  hysteresis  was  observed  in  the  magnetic  transitions 
of  Gd  and  Sm. 

I.  INTRODUCTION 

Magnetic  susceptibility  measurements  were  made  on  several 

polycrystalline  rare  earth  metals  of  99.  8%  purity  using  the  low  field 

(1) 

a.  c.  mutual  inductance  technique  described  by  Hartshorn  with  the 

(2)  (3) 

bridge  and  cryostat  system  described  by  Jennings  and  Gerstein 

respectively.  Interest  in  the  investigation  of  thermal  hysteresis  in 

these  metals  initially  arose  from  the  discrepancy  between  heat  capacity 

'4) 

work  in  this  laboratory  on  Sm  metal  and  that  of  Roberts.  '  *  Roberts 
reported  a  maximum  in  her  heat  capacity  versus  temperature  curve  at 
13.  8°K.  The  work  in  this  laboratory  showed  a  still  rising  Cp  versus 
temperature  curve  at  IE.  7°K,  the  lower  limit  of  these  measurements. 
Thermal  hysteresis  of  the  magnetic  ordering  process  occurring  in  the 
region  of  the  heat  ca.pacity  maximum  was  considered  to  be  a  possible 
explanation  of  this  discrepancy. 
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In  measuring  the  magnetization  of  solid  materials  using  a.  c. 
methods,  one  must  consider  both  dispersion  and  absorption.  The  dis¬ 
persion,  X'm»  is  associated  with  the  change  in  mutual  inductance  of  the 
sample  coils  per  mole  of  sample  present,  and  the  absorption,  Xum>  is 
related  to  the  a.  c.  losses  associated  with  the  specimen.  The  a.  c. 
measurement  enables  one  to  correlate  the  resistive  changes  of  the  bridg 
with  these  losses.  Experimentally  we  have  observed  that  X"mis  zero 
for  T>T^,  the  antiferromagnetic  ordering  temperature,  from  which  we 
infer  that  eddy  current  losses  are  negligible  in  all  samples  measured 
except  samarium.  The  eddy  current  losses  in  samarium  are  believed 
to  be  due  to  the  large  sample  size.  The  a.  c.  losses  at  T < T^  are 
therefore  primarily  associated  with  magnetic  losses.  Thus  at  T  <T^ 
we  are  looking 'at  an  "effective"  a.  c.  susceptibility  and  can  merely  cor¬ 
relate  changes  in  this  function  with  heat  capacity  and  static  susceptibil¬ 
ity  work.  At  T>Tjsj,  X'm  rapidly  approaches  zero,  magnetic  losses 
become  negligible  and  there  is  reasonable  agreement  between  our 
measured  x'm  and  static  susceptibility  results. 

The  measurements  were  made  at  fields  of  10  ±  1  oersted.  The 

bulk  magnetization  below  the  Curie  point  (T^)  is  attributed  to  reversible 

domain  wall  motion,  such  that  domains  with  components  in  the  direc¬ 

ts) 

tion  of  the  external  field  are  enlarged  and  the  others  are  reduced. 

The  sample  coil  of  the  mutual  inductance  bridge  was  calibrated  using  a 
known  paramagnetic  sait,  Mn(NH^) ^(SO^)^'  SH^O.  This  measurement 
yielded  a  numerical  coil  constant  which  permitted  the  conversion  of  re¬ 
sistive  and  inductive  changes  to  X!'m  and  X'm  respectively.  With  the 
exception  of  Sm,  the  average  sample  size  measured  was  approximately 
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2  x  10" ^  moles.  The  sensitivity  of  the  susceptibility  measurement  with 
this  size  sample  was  such  that  a  value  of  X'm  =  10  emu/mole  could 
be  detected  to  within  one  percent. 

II.  RESULTS  AND  DISCUSSION 

1)  Terbium:  Koehler  et  al.  '  performed  neutron  diffraction 

studies  on  single-crystals  and  polycrystalline  Tb.  They  observed  a 

transition  to  an  antiferromagnetic  state  at  230° K  and  a  subsequent 

transition  to  a  ferromagnetic  state  at  220° K.  They  also  found  that  in 

the  antiferromagnetic  region  the  magnetic  structure  of  Tb  is  helical 

and  the  transition  from  the  a.ntiferromagnetic  to  the  ferromagnetic  state 

(7) 

is  centered  at  224°  K.  Jennings  et  al.  7  reported- a  sharp  anomaly  in 
the  specific  heat  versus  temperature  curve  at  227.  7°K  and  an  anomolous 
region  near  220°  K  which  was  found  to  depend  on  the  thermal  history  of 
the  sample  (Fig.  1). 

The  present  results  are  shown  in  Fig.  1.  The  measurements  in¬ 
dicated  a  transition  to  an  ordered  magnetic  structure  at  230.  5°K,  in 
agreement  with  other  workers.  In  the  region  of  ferromagnetic  order 
we  observed  a  four  degree  thermal  hysteresis  effect  upon  heating  and 
cooling.  Data  taken  while  heating  the  sample  showed  a  maximum  in 
X*  at  219-  9°K  and  data  taken  while  cooling  showed  a  maximum  at 
215.  7°K.  The  maximum  on  cooling  always  occurred  at  the  same  tem¬ 
perature  regardless  of  the  time  elapsed  on  cooling.  The  heating  curve 
was  reproduced  only  if  the  sample  was  cooled  to  77° K  and  allowed  to 
stand  for  several  hours.  Otherwise,  various  curves  would  be  obtained, 
ranging  from  the  original  heating  curve  to  the  curve  representing  the 
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sample  being  cooled  from  room  temperature.  Above  the  Neel  point  at 
230.  5°K,  the  heating  and  cooling  curves  coincided.  The  thermal 
hysteresis  effect  was  also  observed  in  the  Xum  versus  temperature 
curve  in  the  region  of  ferromagnetic  ordering.  The  minimum  in  the 
X'm  curve  at  224° K  occurs  in  the  small  but  finite  temperature  region 
where  Koeler  et  al.  ^  report  that  the  interlayer  turn  angle  remains 
constant. 

An  extrapolation  of  the  l/x'm  versus  T°K  curve  above  230°K  to 

l/^m  =  ®  yielded  a  paramagnetic  Curie  point  of  235.  6°K.  The  a-,  b- 

axis  crystals  of  Heglund  et  al.  '  yielded  a  paramagnetic  Curie  point  of 

239°K.  From  their  plot  of  l/x  versus.  T°K  for  a-,  b-axis  crystals  a 

-3 

molar  susceptibility  of  185  x  10  emu.  was  calculated  at  300° K.  Our 
results  yielded  a  molar  susceptibility  at  300° K  of  166  x  10  emu. 

2)  Dysprosium:  The  results  of  the  present  work  are  shown  in 
Fig.  2.  Also  shown  are  the  results  of  heat  capacity  measurements  on 
Dy  metal  performed  in  this  laboratory  by  Griffel  et  al.  K  1  In  this  heat 
capacity  work  two  peaks  are  evident,  the  higher  of  which,  at  174°K, 
was  termed  reproducible  and  did  not  exhibit  thermal  hysteresis.  The 
lower  peak  at  83.  5°K  does  exhibit  a  dependence  on  the  thermal  history 
of  the  sample.  Our  susceptibility  measurements  show  a  transition  at 
179°K,  the  temperature  at  which  other  workers  find  the  onset  of  anti¬ 
ferromagnetic  order.  In  the  region  of  ferromagnetic  order  we  observe 
a  thirteen  degree  thermal  hysteresis  effect.  Data  taken  while  cooling 
the  sample  from  room  temperature  exhibited  a  maximum  in  the  suscepti¬ 
bility  versus  temperature  curve  at  70.  5°K.  This  cooling  maximum  was 
reproducible  regardless  of  what  temperature  above  the  ferromagnetic 
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ordering  temperature  cooling  was  started.  Data  taken  while  heating 
the  sample  from  4.  2°K  exhibited  a  maximum  in  the  susceptibility 
versus  temperature  curve  at  83° K.  This  curve  was  reproducible  pro¬ 
viding  the  sample  was  allowed  to  remain  at  4.  2°K  for  approximately 
one  hour  before  heating.  Above  the  Neel  point  at  179°K,  the  heating 
and  cooling  curves  coincided.  The  susceptibility  data  above  179°K 
yielded  a  paramagnetic  Curie  point  of  152.  5°K  which  agrees  closely 
with  the  results  of  Behrendt  et  al.  and  Arajs  and  Colvin.  Our 
results  also  yielded  a  molar  susceptibility  at  300° K  of  115  x  10  emu. 

The  molar  susceptibility  calculated  by  Arajs  and  Colvin  at  300° K  was 
_3 

1  00  x  10  emu. 

3)  Holmium:  The  results  of  this  work  and  the  heat  capacity  work 

(12) 

of  Gerstein  et  al.  are  shown  in  Fig.  3.  Ho  exhibits  two  peaks  in 
the  heat  capacity  curve.  The  higher  one  at  131.  6°  K  is  reproducible 
and  exhibits  no  dependence  on  the  thermal  history  of  the  sample.  The 
lower  peak  at  19.  4° K  does  exhibit  thermal  hysteresis,  in  that  the  height 
of  the  peak  varies  depending  on  the  thermal  history  of  the  sample.  The 
interval  which  exhibited  thermal  hysteresis  ranged  from  17  to  21° K. 

The  magnetic  susceptibility  measurements  reported  here  also  exhibited 
thermal  hysteresis  in  the  region  of  ferromagnetic  order  near  20° K. 

As  shown  in  Fig.  3,  the  transition  from  the  paramagnetic  to  the  anti¬ 
ferromagnetic  state  was  observed  to  occur  at  129.  3°K.  From  about 
70°  below  and  at  all  temperatures  above  this  point,  the  heating  and 
cooling  curves  coincided.  The  antiferromagnetic -ferromagnetic  transi¬ 
tion,  however,  exhibited  a  2.  7°  thermal  hysteresis  effect  between  heat¬ 
ing  and  cooling  curves.  Before  heating,  the  sample  was  allowed  to 
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remain  at  4.  2°K  for  at  least  one  hour,  thereby  making  the  heating 

curves  reproducible.  The  cooling  curves,  like  those  for  Dy  and  Tb, 

were  reproducible  regardless  of  the  time  elapsed  on  cooling.  The 

maximum  in  the  cooling  curve  occurred  at  18°K  and  the  maximum  upon 

heating  occurred  at  20.  7°K.  The  plots  of  X"m  versus  T°K  exhibited 

maximums  at  very  nearly  the  same  temperature  as  the  corresponding 

X1  plot.  X"  also  approached  zero  much  more  rapidly  in  Ho  than  in 
0  * 

the  other  metals  studied. 

An  extrapolation  of  the  l/x'm  data  to  zero  at  temperatures  con¬ 
siderably  above  the  Neel  point  yielded  a  paramagnetic  Curie  point  of 

/i  o  \ 

86° K.  Strandberg  et  al.  1  calculated  the  magnetic  susceptibility  of 

polycrystalline  Ho  from  susceptibility  measurements  on  a-  and  c-axis 

crystals.  The  molar  susceptibility  at  200° K  calculated  from  their  data 
-3 

is  127  x  10  emu.  Our  a.  c.  measurements  yielded  a  molar  susceptibil¬ 
ity  at  200°  K  of  138  x  10  ^  emu. 

4)  Thulium :  Previous  work  on  the  magnetic  and  thermal  proper¬ 
ties  of  Tm  indicate  that  the  metal's  behavior  is  out  of  the  ordinary 
with  respect  to  the  preceeding  metals  investigated  in  this  work,  namely 
Tb,  Dy  and  Ho.  Neutron  diffraction  studies  by  Koehler  et  al. 
report  that  polycrystalline  Tm  undergoes  a  paramagnetic  to  antifer¬ 
romagnetic  transition  at  53° K  and  a  subsequent  transition  to  a  "fer- 
rimagnetic"  type  of  antiphase  domain  structure  at  38° K.  A  simple 
oscillating  z-component  antiferromagnetic  structure  is  said  to  exist 
between  38  and  53° K.  In  the  "ferrimagnetic"  state  the  net  moment  is 
reported  to  be  parallel  to  the  c-axis. 
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The  results  of  this  work  are  shown  in  Fig.  4.  Heat  capacity- 

measurements  by  Jennings  et  al.  are  included  for  compa.rison.  It 

is  interesting  to  note  that  the  heat  capacity  exhibits  no  peak  below  55° K. 

The  plots  of  X'm  versus  T°K  in  Fig.  4(a)  indicate  that  a  paramagnetic 

to  antiferromagnetic  transition  occurs  at  57.  5°K  The  heating  and 

cooling  curves  above  this  temperature  coincide  and  no  thermal  hysteresis 

is  observed  at  the  Neel  point.  This  is  in  contrast  to  the  lambda-type 

anomaly  in  the  heat  capacity  at  55°  K  which  was  reported  to  exhibit 

thermal  hysteresis.  A  2.8°  thermal  hysteresis  is  observed,  however, 

in  the  region  of  the  antiphase  domain  structure.  The  cooling  curve 

reaches  a  maximum  in  x‘m  at  27.  8°K  and  the  heating  curve  reaches  a 

maximum  in  X’m  at  30.  6°K.  As  in  the  previous  metals  discussed,  the 

heating  and  cooling  curves  were  reproducible.  Fig.  4(b)  illustrates  the 

thermal  hysteresis  effect  in  plots  of  X"m  versus  T°K.  The  elbow  in 

the  curve  at  37.  5°K  agrees  well  with  the  onset  of  "ferrimagnetism"  as 

(14) 

reported  by  Koehler  et  al  in  their  neutron  diffraction  investigation. 

An  extrapolation  of  l/x'm  zero  yielded  a  paramagnetic  Curie 
point  of  17.  5°K.  This  is  in  close  agreement  with  the  results  of  Arajs^"^ 

/  1  Q\ 

Rhodes  et  al.  '  '  obtained  a  molar  susceptibility  at  100°K  of  approxi- 

_  3 

mately  95  x  10  emu.  Our  results  yielded  a  molar  susceptibility  at 
100°K  of  92  x  10"3  emu. 

Mutual  inductance  magnetic  susceptibility  measurements  were 
also  made  on  Gd  and  Sm.  No  thermal  hysteresis  was  observed  near 
the  ferromagnetic  ordering  temperature  (292.  5°K)  in  Gd.  Likewise  no 
thermal  hysteresis  was  observed  in  the  maximum  obtained  at  14.  8°K 
in  Sm.  It  has  been  suggested  that  crystallographic  distortion  near  the 
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Curie  point  in  the  metals  exhibiting  thermal  hysteresis  may  account 

(1 91 

for  these  effects.  '  The  extraordinarily  large  magnetostriction 
values  reported  for  Dy  tends  to  correlate  with  the  large  thermal 
hysteresis  of  this  metal  observed  in  this  work. 
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Fig.  2.  X'  •  X"m  and  C  versus  T*K  for  Dy  metal. 


Fig.  3.  X*  *  X"m  and  C  versus  T'K  for  Ho  metal. 
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Fie.  4.  a)  X'  and  C  versus  T°K  for  Tm  metal. 
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b)  X"  versus  T°K  for  Tm  metal. 
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DISTRIBUTION  OF  MAGNETIC  MOMENT  DENSITY  IN  Tb* 
f 

0.  Steinsvoll  ,  G.  Shirane,  R.  Nathans 
Brookhaven  National  Laboratory,  Upton,  New  York 
H.  A.  Alperin  and  S.  J.  Pickart 

U.  S.  Naval  Ordnance  Laboratory,  Silver  Spring,  Maryland  and 
Brookhaven  National  Laboratory,  Upton,  New  York 

Abstract 

The  magnetic  form  factor  of  metallic  terbium  has  been  measured  at 
4.2°K  by  scattering  of  polarized  neutrons  with  wavelength  X  =  0.8A 
from  single  crystal  samples.  Using  a  split  coil  superconducting  mag¬ 
net  giving  a  field  of  40K  Oe  across  the  crystal  it  was  possible  to 
magnetize  the  samples  close  to  saturation  along  both  the  a-  and  the  b- 
directions,  the  latter  being  the  direction  of  easy  magnetization.  The 
over-all  agreement  with  theoretical  calculations  for  the  Tb3*  ion 
form  factor  is  satisfactory.  The  experimental  points  were  found  to  lie 
on  a  smooth  curve  which  is  consistent  with  the  theoretical  prediction 
that  the  terbium  magnetic  moment  density  has  a  disc-shaped  distribu¬ 
tion  normal  to  the  resultant  moment  direction. 

INTRODUCTION 

The  determination  of  the  neutron  magnetic  form  factors  of  the 
magnetic  elements  is  of  great  value  because  their  Fourier  transforms 
will  give  information  about  the  spatial  distribution  of  the  magnetic 
moment  densities  in  the  atoms.  Much  experimental  and  theoretical  work 
in  recent  years  has  therefore  been  concerned  with  the  determination  of 
the  magnetic  form  factor  for  ions  of  the  transition  elements.  The 
form  factors  for  the  3d  metals  have  been  measured  with  polarized 
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,  neutrons  with  considerable  accuracy.  However,  to  date  no  comparable 

i  measurements  have  been  performed  on  the  rare  earth  metals .  The  diffi¬ 

culty  in  making  such  measurements  is  twofold.  First,  good  single  cry¬ 
stal  samples  must  be  available  and  secondly,  a  magnetic  field  strong 
enough  to  obtain  saturation  must  be  provided  without  depolarizing  the 
incident  neutron  beam. 

An  excellent  single  crystal  of  the  element  terbium  was  put  at  our 
disposal  by  Prof.  F.  H.  Spedding,  Prof.  S.  Legvold,  and  Dr.  L.  Sill  of 
the  Ames  Laboratory.  A  special  superconducting  magnet  was  constructed 
which  provided  a  field  of  40K  Oe  across  the  sample.  The  problem  of 
depolarization  of  the  neutron  beam  through  this  magnet  was  avoided  by 
a  method  to  be  described  below. 

Metallic  terbium  crystallizes  with  a  h.c.p.  structure.  Below 
220°K  it  has  an  ordered  ferromagnetic  state  with  the  b-axis  the  direc¬ 
tion  of  easy  magnetization.  A  moment  of  9.34  at  4.2°K  was  found  by 
D.  E.  Hegland  et  al^  from  magnetization  measurements. 

EXPERIMENTAL 

The  polarized  neutron  technique  has  been  described  in  detail  else- 
2 

where.  Briefly,  the  experiment  consists  of  measuring  the  ratio  be¬ 
tween  two  scattered  intensities  in  the  Bragg  peaks  from  the  sample 
when  the  neutrons  are  polarized  parallel  or  antiparallel  to  the  mag¬ 
netic  field.  This  so-called  flipping  ratio  R  is  related  to  the  ratio 
Y  between  the  magnetic  and  nuclear  scattering  amplitudes,  p  and  b  res¬ 
pectively,  in  the  following  manner: 

r  =  <d 

The  magnetic  scattering  amplitude  is  proportional  to  the  magnetic  form 
factor  of  the  atoms  in  the  sample.  Several  corrections  have  to  be 
applied  to  formula  (1)  in  an  actual  experiment. 

I 

The  original  single  crystal  of  terbium  was  cut  into  smaller 
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pieces  which  could  be  magnetized  along  different  crystallographic 
directions.  In  addition,  the  small-sized  crystals  were  needed  in 
order  to  avoid  secondary  extinction  effects.  This  will  be  discussed 
below.  One  piece  3mm  x  3mm  x  9mm  in  size  was  cut  with  its  long  di¬ 
mension  along  the  b-axis.  Another  smaller  piece  (1mm  x  1.4mm  x  4mm) 
had  its  long  dimension  in  the  a-direction.  The  measurements  were 
performed  at  4.2°K  with  the  crystals  magnetized  along  their  long  di¬ 
mensions  to  minimize  demagnetization.  A  neutron  wavelength  of  0.8A 
was  used  which  permitted  measurements  of  the  form  factor  up  to  the 
value  ^  =  1.1A"1. 

A  schematic  view  of  the  split  coil  superconducting  magnet  is 
shown  in  Fig.  1.  A  doughnut- shaped  region  in  the  space  between  the 
coils  has  a  very  strong  depolarizing  effect  upon  the  neutron  beam. 
This  comes  about  because  the  magnetic  field  changes  direction  in  this 
region,  and  over  a  short  distance  all  magnetic  components  effectively 
go  to  zero.  This  depolarizing  effect  could  be  seen  by  studying  the 
flipping  ratio  at  different  heights  in  the  reflected  beam  from  an 
analyzing  crystal  placed  behind  the  magnet.  The  two  coils  have 
separate  power  supplies  in  order  that  the  current  in  the  coils  can  be 
adjusted  independently.  By  letting  the  lower  coil  carry  less  current 
than  the  upper  one  (see  figure  lb)  this  depolarizing  region  can  be 
forced  to  move  downwards  and  out  of  the  collimated  incident  neutron 
beam.  The  polarization  vector  of  the  primary  neutron  beam  still  has 
to  perform  an  adiabatic  turn  with  the  magnetic  field  in  the  region 
where  the  field  changes  direction  between  the  coils  above  the  bad 
spot.  Here,  however,  a  horizontal  field  component  guides  the  neutron 
polarization  across  the  turning  region,  and  the  polarization  performs 
an  almost  perfect  turn  with  the  magnetic  field.  A  polarization  of 
98%  could  be  measured  in  the  beam,  after  it  had  passed  through 
the  magnet.  In  fact,  the  polarization  had  then  been  turned  twice. 

If  the  sample  is  not  completely  saturated  additional  depolari¬ 
zation  will  occur  in  the  crystal  itself.  A  basic  requirement  in  the 
experiment  is  that  the  sample  be  magnetized  normal  to  the  scattering 
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Figs,  la  and  lb:  1  and  2:  Superconducting  coils.  3:  Crystal  sample. 
4:  Incident  neutron  beam.  5:  Scattered  neutron  beam.  6:  Graph 
showing  flipping  ratio  R  when  measured  as  a  function  of  height  across 
the  neutron  beam.  7:  Polarization  vector  P  of  the  neutron  beam  (P 
polarization  vector  when  the  R.  F.  flipper  is  on).  8:  Return  field 
from  the  superconducting  magnet.  9:  Depolarizing  region.  Fig.  la 
shows  the  situation  when  both  coils  carry  the  same  current.  The  de¬ 
polarizing  region  is  situated  in  the  middle  between  the  coils.  In 
Fig.  lb  the  two  coils  carry  different  currents,  thus  the  depolarizing 
region  is  forced  downwards. 
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vector.  The  strong  field  of  the  superconducting  magnet  made  it  pos¬ 
sible  to  achieve  a  magnetization  close  to  saturation  along  the  a-  and 
b-  directions;  the  small  depolarization  still  present  reduced  the  effec 
tive  beam  polarization  in  the  sample  to  93%. 

Intense  reflections  arising  from  the  large  magnetic  scattering 
-12 

amplitude  (2.51  x  10  cm  in  the  forward  direction)  made  it  necessary 
to  use  two  crystals  of  different  sizes  in  order  to  eliminate  the  ef¬ 
fect  of  secondary  extinction.  A  small  crystal  (the  a-axis  crystal) 
was  used  to  study  the  predominantly  strong  low  angle  ref lections, while 
a  larger  b-axis  crystal  was  used  for  the  weaker  high  angle  reflections. 
Integrated  intensities  measured  on  the  a-axis  crystal  when  plotted 
against  calculated  intensities  showed  no  systematic  trend  indicative 
of  secondary  extinction.  Since  the  (006)  reflection  measured  on  the 
b-axis  crystal  yielded  a  value  of  the  flipping  ratio  in  agreement  with 
the  same  reflection  measured  on  the  extinct ion- free  a-crystal  it  is 
therefore  safe  to  conclude  that  the  even  weaker  high  angle  reflections 
from  the  b-crystal  must  be  free  of  extinction.  Confidence  in  the  final 
form  factor  measurements  is  sustained  by  noting  that  the  strong  reflec¬ 
tions  ((002),  (004),  (006),  etc.)  and  the  weak  reflections  ((100), 
(102),  (200),  etc.  )  all  fall  on  the  same  smooth  curve.  The  presence 
of  simultaneous  Bragg  reflections  was  minimized  by  using  long  pillar¬ 
shaped  samples. 

RESULTS 

The  normalized  magnetic  form  factor  values  were  obtained  using  a 

-12 

value  of  the  nuclear  scattering  amplitude  of  b  =  0.76  x  10  cm  and  a 

value  of  S.34  Ug  for  the  moment  on  the  Tb  atoms.  The  results  are 

given  in  Table  I  and  Figure  2.  For  comparison  the  magnetic  form 

3  3+ 

factor  computed  by  M.  Blume  for  the  Tb  ion  in  an  external  magnetic 
field  has  also  been  drawn. 

A  comparison  of  the  theoretical  value  for  the  saturation  moment 
(S  =  3,  L  =  3,  so  2S+L  =  9  Bohr  magnetons)  with  the  measured  value 
leads  one  to  conclude  that  the  exchange  energy  is  probably  large 
compared  to 
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TABLE 


hki 

(T9)*-1 

f 

hkf  1 

f  sine^ 
k  X  J 

A'1 

f 

100 

.160 

.97 

± 

.02 

224 

.657 

.183 

± 

.002 

002 

.175 

.86 

± 

.02 

008  ' 

.702 

.148 

± 

.001 

101 

.182 

.84 

± 

.02 

118 

.755 

.119 

± 

.002 

102 

.237 

.74 

± 

.03 

226 

.765 

.115 

± 

.001 

103 

.308 

.59 

± 

.03 

330 

.833 

.080 

+ 

.001 

200 

.320 

.65 

± 

.02 

332 

.851 

.075 

± 

.001 

201 

.332 

.61 

± 

.01 

0010 

.878 

.066 

± 

.001 

004 

.351 

.56 

± 

.01 

228 

.895 

.059 

± 

.001 

203 

.415 

.46 

± 

.02 

334 

.904 

.056 

± 

.002 

105 

.467 

.39 

± 

.02 

1110 

.921 

.053 

± 

.005 

300 

.481 

.40 

± 

.03 

336 

.985 

.035 

± 

.001 

302 

.512 

.36 

± 

.03 

2210 

1.04 

.022 

± 

.001 

006 

.527 

.311  ±.006 

1112 

1.09 

.016 

dt 

.001 

222 

.582 

.261  ±.005 

338 

1.09 

.014 

± 

.001 

116 

.596 

.239  i 

r.004 

440 

1.11 

.012 

± 

.002 

The  values  of  the  normalized  magnetic  form  factor  f  have  been 
computed  from  the  experimental  data  assuming  an  effective  polari¬ 
zation  of  93%  and  a  flipping  efficiency  of  997«. 


the  crystal  field  energy  and  thus  Blume's  calculation  assumes  no 
mixing  of  the  Mj  =  6  ground  state  with  higher  states.  This  results 
in  a  disc-shaped  spatial  distribution  of  the  magnetic  moment  density 
which  can  be  oriented  by  the  external  field,  and  justifies  the 
plotting  of  the  a-axis  and  b~axis  data  on  one  curve.  The  geometry  of 
our  experimental  arrangement  (scattering  vector  always  perpendicular 
to  the  direction  of  the  applied  field)  however,  ensures  that  only  the 
rotationally  symmetric  projection  of  the  disc  will  be  observed.  This 
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Fig.  2.  The  normalized  form  factor  of  metallic  terbium  as  a  function 


The  smooth  curve  is  due  to  a  calculation  of  M. 


Blume  as  described  in  the  text. 


prediction  is  confirmed  by  the  fact  that  all  experimental  points  for 
large  values  of  sin9/X  lie  on  a  smooth  curve. 

The  experimental  points  are  somewhat  lower  than  the  theoretical 
curve  at  high  angles.  This  same  tendency  has  also  been  noted  for 
other  rare  fearth  ions  when  their  magnetic  form  factors  have  been 
measured  and  compared  with  appropriate  theoretical  calculations. 


r 
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These  measurements  for  terbium  metal  lie  somewhat  higher  than  the 
values  which  can  be  obtained  from  paramagnetic  scattering  (to  sinP/X  = 
0.6)  by  A'coji^  on  TbC^ . 

The  data  were  not  accurate  enough  at  small  scattering  angles  to 
see  any  possible  scattering  effects  due  to  polarized  conduction  elec¬ 
trons  , 
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MAGNETOSTRICTION  IN  THE  HEAVY  RARE  EARTH  METALS 

B.  F„  De Savage  and  A.  E.  Clark 

U.  S.  Naval  Ordnance  Laboratory, White  Oak 
Silver  Spring,  Maryland 

ABSTRACT 

We  have  found  that  the  magnetostrictive  strains  in 
single  crystals  of  the  heavy  rare  earth  metals  Tb,  Ly,  and 
Er  are  extraordinarily  large  (10“^  to  10”  )  at  low  temp¬ 
eratures.  When  these  strains  are  combined  group  theo¬ 
retically  with  spin-operators  which  transform  according  to 
the  irreducible  representations  of  the  hexagonal  group,  it 
is  observed  that  the  lowest  order  magnetostrictive  strains 
arise  primarily  from  single-ion  terms  in  the  interaction 
Hamiltonian  and  hence  have  a  temperature  and  field  de¬ 
pendence  given  by  a  simple  function  of  the  magnetic  moment: 

X(T)  -  X(0)  ^/2  {f"1[m(T,H)} 

Here  X  [m(T,H)]  is  the  inverse  langevin  function  of  the 

A 

moment  and  I,  ,  .  is  a  red.uced  hyperbolic  Bessel  function 

l  +  a* 

1  8 

of  order  i  +  ^  *  *here  t  is  the  degree  of  the  spin  oper¬ 
ators.  Assuming  the  deformable-ion  model,  Tsuya,  Clark, 
and  Bozorth  predict  the  0*K  values  of  these  strains  x(0) 
from  the  strained  charge  distribution  of  the  conduction 
electrons.  We  find  the  extrapolated  0eK  values  of  the 


.1 


shear  magnetostriction,  taking  I}y  as  100,  to  be  117  >  100, 

28,  and  -44  for  Tb,  Dy,  Ho,  and  Er,  respectively,  which 
are  to  be  compared  with  the  theoretical  estimates  of  115> 
100,  36,  and  -53. 

INTRODUCTION  \ 

Callen  and  Callen1  have  calculated  the  temperature  de¬ 
pendence  of  the  magnetostriction  in  crystals  of  hexagonal 
close  packed  symmetry.  Their  procedure  is  to  write  the  Ham¬ 
iltonian  of  the  system  for  small  strains  as: 

H  -  He  +  \  +  (1) 

where  He  is  the  elastic  energy,  H^  is  the  magnetic  energy 
of  the  unstrained  crystal,  and  is  the  magnetoelastic 
energy.  Callen  and  Callen  express  as: 

\e  “  5,J,i  'l'3  l  •  <2> 

In  this  expression  (S)  are  polynomials  of  degree  i  in 

spin  S  which  transform  according  to  the  i  basis  function 
of  the  yth  irreducible  representation,  are  symmetry 

strain  components  and  are  temperature  independent 

magnetoelastic  coupling  coefficients.  The  index  J  enumer¬ 
ates  strains  of  the  same  irreducible  representation. 

Minimizing  the  free  energy  with  respect  to  the  strains 
and  converting  the  symmetry  strains  to  conventional  strains 
.  referred  to  a  cartesian  coordinate  system  coinciding  with 
the  hexagonal  a-,  b-,  and  coaxes,  the  magnetostriction  is 
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given  by:' 


?  1*1  ,8,Bi  ’  ^,(B*  +  »•)(«*  -  ,l)  +  (or*  -  M 

1>J  i>J  i  J  i  x  y  z  9  a  z  z  a 

+  xY’*{i(e!c  ‘  8y>  («4  "  “yl  +  2My«x“y}  (5) 


+  X«>*(bx«x  +  By0(y)  Szoz 


the  or's  and  p's  are  direction  cosines  of  the  magnetization 
direction  and  measurement  direction,  respectively.  The 
coefficients  are  temperature  dependent  coefficients 
which  are  given  by: 

XT  -  V.  r*"1  (■)]•  w 

In  this  expression  ij**  (o)  Is  the  0*K  value  of  the  magneto- 
striction  and  I  /  [l~x  (m)]  is  a  hyperbolic  Bessel  function 
whose  argument  is  the  inverse  Langevin  function  of  the  re¬ 
duced  magnetization.  The  function  I  /  reduces  to  m®  at 

B/  8 

low  temperatures  and  to  a.m*  at  high  temperatures.  The 
u,a  .  . 

Xj  (0)'s  are  made  up  of  a  magnetoelastic  coupling  constant 
B|j^  divided  by  an  elastic  modulus.  The  superscripts  <*,y, 
and  c  denote  the  three  distinct  modes  of  strain  in  a 
hexagonal  system:  <*,  a  fully  symmetric  strain  in  which 
hexagonal  symmetry  is  preserved,  y,  a  shearing  strain  in 
the  basal  plane,  and  **  a  shearing  strain  in  a  plane  con¬ 
taining  the  £~axls.  The  coefficients  subscripted  1  and  2 
represent  a  change  in  the  basal  plane  diameter  and  the 
axis  length,  respectively. 


Equation  (4)  applies  only  to  one-ion  interactions  in 
which  there  is  coupling  between  the  strain  field  and 
single  spins.  A  two-ion  theory  has  been  derived  by  Callen 
and  Callen,1  however,  it  appears  that  for  the  heavy  rare 
earths  a  single-ion  theory  is  sufficient  to  predict  the 
correct  temperature  dependence  of  the  magnetostriction3 
components  of  Eq.  (3).  For  two  lower-order  components  not 
included  in  Eq.  (5),  the  two-ion  terms  are  required. 

Tsuya,  Clark,  and  Bozorth3  have  calculated  the  magni¬ 
tudes  of  the  magnetostriction  in  the  heavy  rare  earth 
metals  Tb,  Dy,  Ho,  Er,  and  Tm  on  the  basis  of  two  mecha¬ 
nisms:  (l),  the  distortion  of  the  crystalline  electric 
field  due  to  strain  induced  change  of  position  of  the 
nearest  neighbor  positive  ions,  and  (2),  the  distortion  of 
the  electric  field  arising  from  the  3train-perturbed 
charge  distribution  of  the  conduction  electrons.  Using  the 
equivalent  operator  method  of  Stevens  they  find: 

Bj*  -  <*J  (J  -  J-)  <r}>  (5) 

In  this  expression  a  is  a  number  calculated  by  Stevens4 
for  each  of  the  rare  earths.  J  is  the  total  angular 
momentum  in  units  of  h,  and  <r^>  is  the  average  of  the 
square  of  the  radius  of  the  4f  electron  shell. 
is  the  model  dependent  part  of  the  expression.  In  the 
first  mechanism,  W^8  is  a  lattice  sum  containing  terms 


27 


of  the  form  P  «  s(2Z*  -  X*  -  Ya )/2R®  and 
s  n  n  n  n  n 

P  »  s(35Z4  -  30Z8  R*  +  5R4)/8R7.  In  the  second  mechanism, 
4  n  n  n  n  n''  n 

the  magnetostriction  is  estimated  using  the  deformable- 
ion  model  assuming  a  spherical  charge  density.  Under 
these  approximations? 

2Wa's*.IwY’»  *  _iWfi,#  ~  fp(rc)r®  dr. 

8  *  r»J  cc 

and  Wa'p  *  0 
1 

where  p(r  )  is  the  charge  density  of  the  conduction 
electrons. 

It  is  seen  that  the  product  ot  J(«J-l/2)  <rj>  is  common 
to  all  atomic  mechanisms  and  since  varies  slightly 

with  the  number  of  4f  electrons,  a  good  estimate  of  the 
relative  magnitude  of  the  magnetoelastic  coupling 
constants  in  the  rare  earths  is  found  by  a  comparison  of 


the  expression: 


aJ  (J  -  1/2)  <rj> 


After  introducing  the  appropriate  elastic  constants,  the 
0*K  magnetoelastic  constants  of  Eq.  (2)  become: 

x“"(o)  -  25(w°"*--w“'*)o  -  2|(w®’*+ 

1  U  1  3  8  33  U  1  38  13 

X®'8(0)  “  2k(W?'8-  LW*’»)(-2C  )4|(Wff'8+i»W«'*)(C  +C  ) 

8  l)  1  3  8  13  '  D'  1  3  8  ''ll  18^ 


x^8(0) 


KWY>8 


,5  K  -  a<r*->  J(J  -  L) 


KW*>* 


X*"  (0)  -  _* _  ;  D  -  C  (C  +C  )  -  2C  * 

4(5.  "  *S  11  18  18 


Estimating  the  elastic  constants  from  the  room  temper- 
ature  sound  velocities  in  dysprosium,®  we  find: 

x“’*(0)  -  55^-  (15. 7W*"-  12.66W0"*)  10-11 

x“’*(°)  *  zJV-  (14. 8W®'*  +  24.8W®’’  )1CT11  (7) 

S  OH-.J.  1  c 

XY'*(0)  *  KWY>*/9.l4  x  10*1 
X*,#(0)  -  KW€,*/9.76  x  10l* 

EXPERIMENT 

We  have  sought  to  test  this  theory  by  measuring  all 
of  the  magnetostriction  coefficients  in  terbium  in  the 
paramagnetic  region,  where  the  magnetocrystalline  aniso¬ 
tropy  is  negligible.  Ibis,  together  with  previou*  measure¬ 
ments  of  the  magnetostriction  in  Er,»  Dy,»  and  Hb*  Bhould 
give  an  indication  of  the  validity  of  the  approximations 
made  in  the  Tsuya-Clark-Bozorth  theory. 

We  have  used  the  conventional  strain  gage  method  to 
measure  the  magnetostriction  in  terbium.  In  order  to 
separate  the  four  coefficients  experimentally,  it  was 
necessary  to  measure  the  strains  along  the  a-  and  coaxes 
and  a  direction  45*  between  the  two  while  the  magnetic 
field  was  rotated  from  the  measurement  direction  to  a 

direction  perpendicular,  For  example,  x0>*  is  determined 

a 

by  measuring  the  change  in  length  along  the  c-axis  while 
rotating  the  field  from  the  £-  to  the  a-axis.  The  0*K 
magnitude  is  found  by  extrapolating  the  experimental  data 


to  0*K  according  to  the  theoretical  temperature  dependence 
of  Eq.  (4),  This  procedure  requires  that  our  measurements 
be  consistant  with  the  dictates  of  the  Callen  and  Callen 
theory,  namely,  that  at  each  temperature,  the  magnetic 
moment  be  held  constant  during  the  measurement.  In  the 
paramagnetic  region,  the  magnetic  susceptibility  along 
the  £-axia  is  different  from  that  along  the  a-axis.  In 
order  to  maintain  constant  moment  during  measurement,  the 
magnetic  field  must  be  changed  during  rotation  if  the 
initial  and  final  directions  have  different  susceptibilities 

as  for  example  in  the  measurements  of  \a,a  and  \a>*  . 

1  * 

Figure  1  illustrates  this  point.  Ihe  top  curve  is  the 
temperature  dependence  of  x®,#  measured  at  constant  moment. 

9 

The  experimental  points  agree  well  with  the  theoretical 
temperature  dependence.  The  lower  curve  shows  the  magneto¬ 
striction  measured  at  constant  field  and  here  Li/l  is  found 
to  have  the  opposite  sign  and  a  different  temperature 
dependence . 

The  temperature  dependences  of  the  four  coefficients 
are  shown  in  Fig.  2,  which  is  a  log-log  plot  of  the  x's 
versus  the  function  I  /  [x"1 (")] .  Best  agreement  with  the 

Callen  and  Callen  theory  is  found  with  lines  drawn  through 
the  data  points  having  slopes  of  unity  and  intercepts  at 
T  =  0*K  (I8/a  »  1)  of  Xa,a(Q)  =  -  0.26  x  10“*, 

\a’*(0)  =  0.9  X  10“*,  xV'*(0)  =  1.0  x  10~a,  and 
a 
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X6,2(0)  =  1.15  X  10"8. 

DISCUSSION 

Table  I  compares  the  predicted  and  measured  relative 
magnitudes  of  the  magnetostriction  in  the  heavy  rare 
earths  by  making  use  of  Eq.  (5).  The  values  of  <r^>  are 
taken  from  Freeman  and  Watson.7  The  magnetostriction  x 
is  related  to  the  magnetoelastic  coupling  coefficients 
B  by  \  c  -  BC”1 ,  where  C  is  an  elastic  constant.  The 
elastic  constants  have  been  measured  only  for  Dy  and  in 
lieu  of  this,  the  ratio  of  the  polycrystalline  shear 
moduli  have  been  used  to  calculate  the  magnetostriction. 


Table  I 


Element 

a  <r|>a.u. 

xth 

Aexp 

polycrystal 

X 

Tb 

-1/99  .756 

115 

117 

+ 

Vy 

-2/315  .726 

100 

100 

+ 

Ho 

-1/450  .696 

56 

28a# 

+ 

Er 

4/1575  .666 

-55 

-44 

-• 

Tra 

1/99 .638 

-74 

— 

- 

a.  Legvold,  Alstad,  and  Rhyne,  Rhys.  Rev.  Letters  10,  50 
(1965) 

The  magnitude  of  x  as  predicted  by  the  theory  is  given 
relative  to  Dy  taken  arbitrarily  as  100.'  The  experimental 
determinations  of  xY,*(o)  normalized  to  Dy  are  shown  in 
the  table  along  with  the  sign  of  the  magnetostriction  as 
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measured  in  the  polycrystalline  material.  Good  agreement 
between  experiment  and  theory  is  observed  for  the  signs 
and  relative  magnitudes.  The  predicted  sign  reversal  be¬ 
tween  Ho  and  Er  is  observed  in  both  single  crystal  and 
polycrystalline  materials. 

In  Table  II,  our  experimental  results  on  Tb  are  com¬ 
pared  to  those  calculated  on  the  basis  of  the  deformable- 
ion  model  (assuming  a  spherical  charge  density  for  the 
conduction  electrons)  and  those  calculated  from  the  dis¬ 
tortion  of  the  crystalline  field  due  to  the  field.  As  ex¬ 
pected,  because  of  the  screening  of  the  4f  .electrons  by  the 
outer  shells,  the  deforraable-ion  magnetostriction  is  too 
large.  In  the  ®able,  we  have  normalized  the  theoretical 
value  of  xY’8  to  the  experimental  value  of  0.010  by  as¬ 
suming  the  screening  factor  a  to  be  0.75.  This  is  con¬ 
sistent  with  the  value  of  o  *  0.8  estimated  by  Freeman.® 
Excellent  agreement  is  found  between  theory  and  experiment 

for  the  coefficient  \a>* S  the  other  two  coefficients, 

a 

\a*9  and  xe,*>  are  the  correct  sign  and  order  of  raagni- 
i 

tude.  On  the  other  hand,  the  magnetostriction  calculated 
from  the  crystalline  field  is  too  small  and  predicts  the 
wrong  sign  for  x*,a  and  incorrect  ratios  for  x®#*/xY'* 

and  x*'VxY'#-  Hence  it  appears  that  the  crystal  field 

» 

plays  a  minor  role  in  the  magnetostriction  of  the  rare 
earths. and  that  the  dominant  mechanism  is  the  interaction 
between  the  conduction  electrons  and  the  4f  electrons. 

r' 
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Table  II 


X  , 

Experiment  Deformable  Ion 

Crystal  Field 

A1 

-.0026 

-.0045 

-.0017 

x£’* 

.009 

.0088  • 

.0055 

XY’» 

.010 

.010 

.0006 

X*'* 

.0115 

.0094 

-.0025 
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OBSERVATION  OF  OPTICAL  ABSORPTION  BANDS  ASSOCIATED  WITH 


MAGNETIC  ORDER  IN  RARE-EARTH  METALS. 

C.  Chr.  Schiller 


IBM  Zurich  Research  Laboratory 
Riischlikon,  Zurich,  Switzerland 


ABSTRACT 

Spectroscopic  studies  of  Gd-metal  reflexion  methods  have  re¬ 
vealed  an  absorption  band  in  the  infrared  ,  centered  at  0.4  e’V,  which 
appears  in  the  ferromagnetically  ordered  state  below  289°  K.  The 
strength  of  the  absorption  increases  with  increased  order  as  the  tem¬ 
perature  is  lowered. 

A  similar  band  is  present  in  Dy-metal  in  the  spiral-phase  be¬ 
low  178°  K  and  persists  into  the  ferromagnetic  phase  below  85°  K. 

The  nature  of  these  extra  absorptions  is  not  jet  understood.  An 
earlier  explanation  on  the  basis  of  magnetic  superzones  present  in  the 
antiferromagnetic  phases  of  Dy  and  Ho  must  now  be  abandoned, 
since  it  does  not  apply  to  ferromagnetic  order. 

INTRODUCTION 

Optical  spectroscopy  of  metals  yields  valuable  information  on 
the  properties  of  conduction  electrons  and  electronic  energy  bands.  As 
part  of  a  more  extended  investigation  of  the  optical  properties  cf  rare- 
earth  metals,  this  paper  reports  on  experimental  results  which  are  con¬ 
cerned  with  a  special  aspect  of  this  work  :  the  influence  of  maj  netic 
order  on  the  reflectivity  of  Gd,  Dy  and  Ho. 

The  original  purpose  of  the  studies  to  be  reported  was  t>  verify 
theoretical  predictions  by  Miwa  (1962).  He  noticed  that  optical  i  osorp- 
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tion  bands  should  be  observable,  associated  with  periodic  spin  order 
which  is  found  in  many  rare  earth  materials. 

It  had  been  observed  earlier  that  the  electrical  resistivity  of 
Tb,  Dy,  Ho  and  Er  shows  profound  anomalies  in  the  antiferromagnetic 
phases  of  these  metals,  which  are  characterized  by  periodic  spin  ar¬ 
rangements.  In  order  to  explain  these  anomalies,  Mackintosh  (1962)  and 
Miwa  (1962)  suggested  the  following  model  : 

The  periodic  arrangement  of  localized  ion  spins  interacts  with 
the  conduction  electrons  through  the  well-established  s-f  exchange  inter¬ 
action.  Since  the  periodicity  is  different  from  that  of  the  lattice,  addi¬ 
tional  Brillonin  zones  are  formed,  and  new  energy-gaps  are  introduced 
into  the  conduction  band.  The  resulting  change  in  Fermi  surface  is 
reflected  by  a  striking  increase  in  electrical  resistivity. 

Calculations  on  the  basis  of  this  model  by  Miwa  and  in  more 
detail  by  Elliott  and  Wedgwood  (1963,  1964)  are  in  good  agreement 
with  the  experimental  data.  Superzone  boundaries  should,  however,  not 
only  influence  the  electrical  resistivity  but  also  give  rise  to  optical 
absorption  from  electron  transitions  across  the  new  gaps.  The  gaps  are 
of  the  order  Eg  =  VS  M(T)  If  (V)  |  where  V  is  an  effective  exchange 
energy,  S  the  ion  spin,  M(T)  the  temperature-dependent,  relative 
magnetisation  and  If  (d  a  structure  factor.  Inserting  numerical 
values  and  taking  M(T)  not  close  to  zero,  Eg  becomes  a  few  tenths  of 
an  eV,  which  corresponds  to  optical  absorption  in  the  near  infrared. 

Our  first  reflectivity  measurements  on  Ho  did  indeed  reveal  an  absorp¬ 
tion  band  with  a  peak  at  0.  35  eV  for  the  magnetic  spiral  phase  below 
132°K  (Schiller,  1964),  which  seemed  to  fit  the  predictions  rather  well. 

Further  extension  of  this  work  to  Dy  and  Gd,  will  now  be  pre¬ 
sented,  which  calls  for  some  modifications  of  the  original  ideas. 
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EXPERIMENTAL 


All  measurements  were  done  on  vacuum  deposited,  polycrystal¬ 
line  films,  about  3000  &  thick,  which  provided  the  necessary  optical  sur¬ 
faces  in  a  convenient  manner.  In  order  to  obtain  film  samples  with  elec¬ 
trical  and  magnetic  properties  close  to  those  of  bulk  material,  all  depo¬ 
sitions  were  carried  out  under  ultra-high  vacuum  conditions.  This,  to¬ 
gether  with  an  annealing  procedure,  produces  films  of  the  required 
quality,  as  we  have  demonstrated  earlier  (Schiller,  1963).  The  samples 
were  kept  in  the  same  vacuum  throughout  the  measurements  so  that  no 
problem  with  surface  contamination  arose.  A  dewar  finger  to  which  the 
substrates  were  attached  allowed  cooling  of  the  specimen  down  to  78°K. 

Reflectivity  at  nearly  normal  incidence  was  determined  by  com¬ 
paring  with  a  silver  film  deposited  simultaneously  and  close  to  the 
rare-earth  film.  The  vacuum  chamber  was  fitted  with  a  sapphire  window 
so  that  measurements  between  0.  25  and  5  eV  of  photon  energy  were 
possible. 

RESULTS 

First  the  reflectivity  R  of  the  sample  at  room  temperature  was 
found  by  comparing  with  the  reflection  from  the  Ag  film  for  which  very 
accurate  absolute  reflectivity  data  are  available  from  the  work  of 
Bennet  et.  al  (1965).  Afterwards,  reflexion  vs.  temperature  curves 
were  taken  for  a  large  number  of  fixed  photon  energies  between  0.  25 

cincl  9  gV. 

Figure  1  shows  an  R  vs.  T  plot  for  Dy.  One  notices  an  increase 
in  R  with  decreasing  temperature.  This  corresponds  to  a  reduction  in 
phonon  scattering  for  the  coziduction  electrons  which  dominate  the  optical 
properties  in  this  energy  range.  At  179°  K,  the  ordering  temperature 
for  the  magnetic  spiral  structure,  there  is  a  sudden  change  in  slope, 
and  R  now  decreases  as  the  temperature  is  lowered  further,  indicating 


additional  absorption.  Similar  curves  were  found  earlier  for  Ho. 

It  is  rather  surprising,  however,  that  there  is  no  evidence  of  a 

discontinuity  accompanying  the  transition  to  the  f *  romagnetic  state  at 

T  =  85°  K.  This  is  in  contrast  to  the  behaviour  of  the  electrical  resis- 
c 

tivity,  where  a  distinct  drop  at  indicates  that  the  superzone  bound¬ 
aries  and  energy  gaps  vanish. 

An  obvious  choice  for  an  experiment  which  can  clarify  this  point 
is  Gd,  which  is  only  simple  ferromagnetic. 

The  result  of  measurements  identical  to  that  on  Dy  was  rather 
unexpected. 

An  R  vs.  T  curve  is  shown  in  Fig.  2  It  is  essentially  of  the 

same  nature  as  for  Dy  except  for  the  different  ordering  temperature 

of  T  =  289°  K. 
c 

This  seems  to  indicate  that  magnetic  order,  regardless  of  its 
periodicity,  plays  the  dominant  role  in  this  effect. 

More  support  is  lent  to  this  point  by  comparing  the  spectral 
dependence  of  the  anomaly,  which  is  plotted  in  Figs.  3  and  4.  As 
previously  for  Ho,  a  resonance-like  dispersion  is  found  with  a  peak 
in  the  0.3  to  0. 4  eV  range.  It  is  remarkable  that  there  are  only  minor 
differences  between  the  curves  for  the  ferromagnetic  state  of  Gd  and 
the  spiral  state  of  Dy. 

An  important  aspect  of  the  absorption  bands  is  the  temperature 
dependence  of  their  position.  It  can  be  seen  from  Figs.  3  and  4  that 
little  shift  towards  lower  energies  is  observed  for  higher  temperatures 
in  Dy,  none  for  Gd.  Such  a  shift  should  however  occur  for  an  energy 
gap  which  varies  with  M(T).  Only  the  magnitude  of  the  anomaly  in¬ 
creases  with  magnetic  order. 


DISCUSSION 


An  explanation  of  the  observed  effect  on  the  basis  of  magnetic 
superzones  must  obviously  be  abandoned  in  view  of  the  results  on  Gd. 

The  observation  of  an  almost  temperature  independent  position  of  the 
band,  unfortunately  also  renders  unlikely  an  association  with  energy 
gaps  produced  by  band- splitting  into  spin-up  and  spin-down  bands  due 
to  magnetic  order,  since  such  gaps,  too,  are  expected  to  be  proportional 
to  M(T). 

An  even  simpler  explanation  one  might  think  of,  is  an  influence 
of  spin-disorder  scattering  of  the  conduction  electrons  which  has  a 
pronounced  effect  on  the  electrical  resistivity  and  is  almost  equally 
reduced  by  periodic  and  ferromagnetic  order. 

However,  our  observations  show  that,  above  the  ordering  point, 
reduction  in  phonon  scattering  at  lower  temperatures  leads  to  an  increase 
in  reflectivity.  So  we  would  expect  even  more  increase  in  reflectivity 
upon  additional  reduction  of  spin  disorder  scattering  by  magnetic  order. 
Instead,  less  reflectivity  is  found  in  the  ordered  state,  which  indicates 
rather  more  absorption. 

In  addition,  the  resonance  like  behaviour  in  the  dispersion  makes 
spin-disorder  an  even  more  unlikely  cause  of  the  effect. 

We  can  only  conclude  that,  so  far,  we  have  been  unable  to 
suggest  a  plausible  explanation  for  our  observations.  Further  studies 
are  in  progress,  using  single  crystal  films  and  polarized  radiation  in 
order  to  take  the  large  anisotropy  of  transport  phenomena  in  the  rare- 
earth  metals  into  account,  which  we  have  neglected  so  far.  It  is  hoped 
to  obtain  some  clues  from  this  extension  of  the  work. 
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ABSTRACT 

This  note  points  out  (1)  that  there  is  a  definitive 
experiment  for  distinguishing  whether  any  part  of  the 
infrared  absorption  (from  0.25  to  0.85  eV)  for  the  heavy 
rare  earth  metals  is  associated  with  the  spiral  periodicity, 

(2)  that,  as  reported  below,  this  experiment  shows  that  for 
dysprosium  the  infrared  absorption  is  completely  indepen¬ 
dent  of  any  periodicity  associated  with  the  spiral,  and 

(3)  that  in  the  light  of  recent  band  calculations,  there  is 
an  alternative  explanation  for  the  occurrence  of  low  fre¬ 
quency  optical  structure  for  the  heavy  rare  earth  metals. 

Recently  Schuler*'  reported  the  appearance  of  structure 
at  about  0.35  eV  in  the  infrared  reflectivity  for  a  holmium 
film  in  the  temperature  range  where  holmium  has  a  spiral 
spin  arrangement  with  periodicity  along  the  c  axis.  Follow- 

n 

Ing  the  ideas  of  Miwa,  Schuler  suggested  that  the  structure 
observed  ia  due  to  the  optical  absorption  corresponding  to 
energy  gaps  in  the  conduction  bands  at  the  magnetic 
Brillouin  zone  boundaries  associated  with  the  spiral 
periodicity.  It  is  the  purpose  of  the  present  note  to  point 
out  (1)  that  there  is  a  definitive  experiment  for  distin¬ 
guishing  whether  this  structure  or  indeed  any  part  of  the 
infrared  absorption  for  the  heavy  rare  earth  metals  is 
associated  with  the  spiral  periodicity,  (2)  that,  as  re¬ 
ported  below,  this  experiment  shows  that  for  dysprosium. 
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which  also  has  spiral  spin  structure,  the  infrared  absorp¬ 
tion  is  completely  independent  of  any  periodicity  associated 

with  the  spiral,  and  (3)  that  in  the  light  of  recent  band 
3  4 

calculations  *  there  is  an  alternative  explanation  for  the 
occurrence  of  such  low  frequency  structure  for  the  heavy 
rare  earth  metals. 

The  definitive  experiment  for  seeing  whether  the  infra¬ 
red  structure  is  associated  with  the  spiral  periodicity 
rests  on  the  fact  that  the  application  of  a  magnetic  field 
in  the  plane  of  the  spiral  can  be  used  to  change  the  spin 
arrangement  to  one  differing  only  slightly  from  ferromag¬ 
netic  alignment  along  the  field.  This  change  occurs  quite 

5 

sharply  at  a  critical  field  H  .  In  general  H  is  not  very 

c  c 

large.  If  the  infrared  structure  is  due  to  the  extra  spiral 
periodicity,  this  structure  can  be  extinguished  by  applying 
a  field  greater  than  H  in  the  plane  of  the  spiral.  This 
effect  will  occur  very  sharply  as  a  function  of  field. 
'Therefore  one  need  only  examine  the  change  in  reflection  or 
absorption  at  a  given  temperature  in  the  spiral  regime  on 
applying  a  field  H  >  H  .  Dysprosium  is  particularly  suit- 
able  for  such  a  study  since  H  is  less  than  10  kilogauss 
throughout  most  of  the  spiral  te  perature  region.  We  have 
therefore  examined  transmission  through  Dy  films  for  a  range 
of  frequency  between  0.25  eV  and  0.85  eV. 

The  measurements  were  made  in  a  cryostat  with  calcium 
fluoride  windows  and  special  pole  pieces  which  allowed  a 
field  of  22  kilogauss  to  be  applied  perpendicular  to  the 
light  path  and  in  the  plane  of  the  film.  The  films  studied 


had  transmissions  in  the  region  of  interest  of  a  few  per¬ 
cent,  and  were  a  few  thousand  Angstroms  thick.  They  were 
prepared  by  evaporation  on  sapphire  substrates  in  a  separate 
vacuum  system.  For  protection  during  the  transfer  to  the 
cryostat  the  films  were  given  a  protective  evaporated  layer 

o 

of  silicon  monoxide  2000  to  3000  A  thick,  which  was  present 
during  the  measurements.  This  would  not  cause  any  change 
in  the  measured  absorption  since  silicon  monoxide  is  trans¬ 
parent  in  the  range  of  the  present  measurements.  Electrical 
measurements  on  auxiliary  samples  indicated  that  such  a  coat 
was  effective  in  preserving  the  metallic  nature  of  the 
dysprosium  film.  No  change  was  observed  in  the  film  proper¬ 
ties  in  several  days  in  the  cryostat.  However,  after  one  to 
two  weeks’  exposure  in  air  the  transmission  increased  sev¬ 
eral  times.  As  a  final  check  on  the  experimental  technique, 
the  difference  in  transmission  between  a  thick  film  and  a 
thin  film  was  studied  and  this  showed  the  same  absorption  as 
indicated  in  the  figure,  A  triangular  five-point  smooth  was 
used  on  the  raw  data. 

The  results  obtained  on  cooling  below  the  Neel  tempera¬ 
ture,  T^  =  179 °K,  are  shown  in  Fig.  1.  The  figure  shows  the 
transmission  at  a  given  temperature  normalized  to  that  at 
182°K  (just  above  T^).  There  is  a  distinct  dip  in  the  trans¬ 
mission  for  the  87 °K  curve  at  about  0.44  e'v  and  at  a  some¬ 
what  lower  frequency  for  123  °K.  (The  lowest  temperature 
shown,  87°K,  is  quite  close  to  that,  *  85°K,  at  which  bulk 
dysprosium  has  a  thermal  transition  to  a  ferromagnetic 
phase.)  Such  structure  is  clearly  absent  for  the  235 °K 
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curve  where  Dy  has  no  magnetic  ordering.  (Structure  is  also 
absent  in  the  corresponding  room  temperature  curve  which 
lies  between  1.3  and  1.5  on  the  scale  of  Fig.  1  for  photon 
energies  between  0.3  eV  and  0.8  eV. )  It  is  difficult  to 
clearly  distinguish  whether  there  is  structure  in  the  fre¬ 
quency  range  observed  at  138 °K.  For  the  spiral  mechanism 
suggested  by  Schuler  the  absorption  structure  should  go  to 
zero  frequency  at  T^,  while  for  the  mechanism  based  on  the 
band  structure  of  the  conduction  electrons  discussed  below 
it  will  probably  shift  to  very  low  frequencies  as  tempera¬ 
ture  increases  toward  TN«  As  the  structure  shifts  to  lower 
frequencies  for  temperatures  increasing  toward  Tn>  it  is 
increasingly  obscured  by  the  free  carrier  contribution  to 
the  absorption.  This  probably  accounts  for  the  difficulty 
in  observing  any  well  defined  structure  at  138 °K. 

There  is  no  measurable  change  in  the  curves  shown  in 
Fig.  1  when  a  field  of  22  kilogauas  is  applied  in  the  plane 
of  the  film.  Ideally  the  experiment  would  be  done  with 
single  crystals.  However,  the  use  of  a  polycrystalline 
film  presents  no  great  obstacle.  For  a  polycrystalline 
sample,  any  given  crystallite  will  h#ve  its  c  axis  at  some 
angle  0  to  H.  The  main  effect  of  the  applied  field  is  to 
produce  an  effective  field  Hsin0  in  the  plane  of  the  spiral 

7 

which  can  distort  the  spiral.  Thes  crystallites  with 
0  >  9  are  essentially  ferromagnetic „  while  those  with 
0  <  0  are  spirals.  Here 

9c  =  sin"1  <Hc/H).  (1) 


Then  the  fraction  of  crystallites  having  a  ferromagnetic 
arrangement  for  a  given  field  is 


f  increases  rapidly  for  H>H  .  For  example,  f  =  0.87  for 

H  =  2H  .  For  T  =  123 °K,  for  dysprosium  H  is  about  5  kilo- 
c  c 

gauss.  Thus  for  an  applied  field  of  22  kilogauss,  H/H  >  4, 
giving  f  >  0.97.  Thus  in  Dy,  through  most  of  the  spiral 
region,  one  can  obtain  almost  complete  ferromagnetic  align¬ 
ment  with  a  field  such  as  the  22  kilogauss  used  in  the 
present  experiment.  Then  just  as  for  single  crystals,  one 
should  be  able  to  obtain  almost  complete  extinction  of  any 
optical  structure  associated  with  the  spiral  periodocity  by 
applying  a  field  of  this  magnitude.  The  fact  that  there  is 
no  change  in  the  infrared  transmission  structure  upon  appli¬ 
cation  of  the  field  clearly  shows  that  the  structure  does 
not  depend  on  the  spiral  periodicity  and  the  resulting 
presence  of  magnetic  Brillouin  zone  boundaries. 


There  is  an  alternative  explanation  for  the  occurrence 

of  such  low  frequency  structure  for  the  heavy  rare  earth 

3  4 

metals.  Recent  work  of  Dimmock,  Freeman  and  Watson  *  indi¬ 


cates  that  the  conduction  electron  bands  for  the  heavy  rare 
earth  metals  closely  resemble  those  of  transition  metals 
rather  than  being  free  electron  like.  This  leads  to  a  large 
density  of  states  at  the  Fermi  level.  The  behavior  for  the 
rare  earth  metals  from  Gd  to  Tm  is  expected  to  be  similar 
in  this  regard.  Since  the  situation  for  the  heavy  rare 
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earths  is  closely  analogous  to  that  in  nickel  and  iron  with 

flat  d  like  bands  lying  close  to  the  Fermi  energy,  it  is 

quite  reasonable  to  expect  the  same  sort  of  low  frequency 

8~ 10  10 

optical  structure  as  occurs  for  nickel  and  iron. 

Indeed  the  frequency  at  which  the  structure  occurs  for  Dy 

and  Ho  is  quite  similar  to  that  for  nickel  and  iron.  The 

optical  structure  in  nickel  is  attributed  to  transitions 

from  a  flat  d  band  to  the  s  band  in  the  vicinity  of  the 
8  11-13 

Fermi  energy.  *  If  a  similar  mechanism  is  the  cause 

of  the  low  frequency  optical  structure  in  Dy  and  Ho,  simi¬ 
lar  structure  should  occur  for  all  the  heavy  rare  earths 
including  gadolinium  which  is  an  ordinary  ferromagnet.  One 
would  expect  the  structure  to  appear  from  very  low  frequen- 
cies  on  cooling  below  T„  ns  the  d  band  for  minority  spin 
electrons  in  the  ordered  magnetic  state  is  driven  increas¬ 
ingly  above  the  Fermi  energy  by  the  exchange  with  localized 
f  electrons.  The  structure  would  shift  to  higher  frequen¬ 
cies  until  the  maximum  d  band  exchange  splitting  is  obtained 
with  thermal  saturation  of  the  raagnetization.  We  also  note 
that  there  is  an  indication  of  en  additional  dip  in  the 
relative  transmission  curves  for  87 °K  and  123 °K  at  about 
0.6  eV.  As  in  the  case  of  nickel,  *  additional  structure 
at  such  frequencies  would  be  quite  reasonable  for  the  model 
just  discussed.  To  conclude  we  point  out  that  the  applica¬ 
tion  of  a  magnetic  field  in  the  way  done  in  the  present 
experiment  would  be  expected  to  have  virtually  no  effect  on 
•  the  infrared  structure  for  such  a  model.  Thus  our  present 

observations  on  Dy  are  completely  consistent  with  what  one 

3  4 

expects  by  considering  the  recent  band  calculations  *  for 
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the  heavy  rare  earth  metals  and  the  known  behavior  of 
Ni. 
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TRANSMISSION  (5)  T°K 
TRANSMISSION  <cD  I62°K 


PHOTON  ENERGY  ~  eV 


Fig.  1.  Relative  transmission  for  dysprosium  films  above  and  below 
the  Neel  temperature  (179 °K). 
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ELECTRON-SPINWAVE  SCATTERING  IN  GADOLINIUM 


B.  Ltithi  arid  H.  Rohrer 
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Rttschlikon,  Zurich,  Switzerland 

ABSTRACT 

The  electrical  resistivity  due  to  magnetic  scattering  in  poly¬ 
crystalline  gadolinium  has  been  determined  in  the  temperature  range 
between  4.  2°  K  and  25°  K.  The  part  of  the  resistivity  due  to  scatter 
ing  by  domains  and  magnetic  impurities  is  practically  independent  of 
temperature,  whereas  that  due  to  spinwave  scattering  follows  a 
t3,  1  +  0.  2  (  jaw  amj  jg  considerably  smaller  than  the  resistivity  due 
to  phonon  scattering.  The  temperature  dependence  of  the  latter  varies 
between  T^‘  ^  and  T^‘  ^  . 


Results  are  presented  here  of  an  experiment  which  determines 
directly  that  part  of  the  electrical  resistivity  due  to  magnetic  scatter¬ 
ing  of  conduction  electrons.  By  applying  a  sufficiently  high  magnetic 
field,  the  magnetic  excitations  will  be  frozen-out  and  the  magnetic 
contribution,  p  m  ,  to  the  electrical  resistivity  will  disappear.  Up  till 
now,  a  Grilneisen  function  for  the  phonon  resistivity  was  usually  sub¬ 
tracted  from  the  total  resistivity  in  order  to  obtain  the  magnetic  part. 

It  will  be  shown  that  this  is  incorrect  at  temperatures  below  25°  K  . 

The  experiment  was  performed  on  Gd  because  it  is  a  rare  earth  met¬ 
al  where  the  s-f  exchange  mechanism  stands  on  somewhat  firm  grounds, 
and  because  it  has  a  low  anisotropy  field  of  about  I  kG  at  helium 
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temperature.  Six  polycrystalline  samples  of  different  origin  were 
studied  with  a  pulsed  magnetic  field  technique  in  the  temperature 
range  between  4  and  25°  K.  In  addition,  the  magnetoresistance  wao 
measured  in  a  longitudinal  field  in  order  to  keep  the  demagnetising 
field  and  the  ordinary  magnetoresistance  low.  Supplementary  mag¬ 
netization  measurements  in  this  geometry  do  indeed  show  a  saturation 
already  at  about  1  kG. 


Figure  1  shows  typical  experimental  results.  Here,  the  meas¬ 
ured  magnetoresistance  is  plotted  against  the  applied  field.  The  ini¬ 
tial  decrease  of  the  electrical  resistivity  is  due  to  the  freezing -out 
of  the  magnetic  scattering  processes  and  the  following  increase  comes 
from  ordinary  longitudinal  mngnetoresistance.  At  low  temperatures, 
the  magnetic  resistance  Ap  ^  will  saturate  in  relatively  low  fields 

4G  kG)  so  that  the  normal  magnetoresistance  Ap  can  be  deter- 

1)  n 

mined  easily.  By  application  of  the  Kohler-rule  the  magnetic  re¬ 
sistivity  for  higher  temperatures  can  then  be  obtained,  and  is  shown 
by  the  dotted  lines  in  Figure  1.  It  should  be  mentioned,  that  the  Kohler - 
rule,  unfortunately,  could  not  be  used  to  compare  the  ordinary  mag¬ 
netoresistance  of  different  samples,  because  most  specimens  showed 
a  pronounced  texture. 


At  higher  temperatures,  saturation  is  no  longer  reached  in  our 
fields  and  an  extrapolation  procedure  has  to  be  applied  in  order  to 
determine  the  total  magnetic  resistivity.  For  high  fields,  a  simple 
exponential  dependence  with  the  spinwave  energy-gap  is  expected 
because  normally  only  one -spinwave  processes  will  be  important.  Such 
an  exponential  plot  of  Ap  does  indeed  give  straight  lines,  and  the 
total  magnetic  resistivity  p^  can  then  be  determined  immediately. 

An  additional  weak  field -and  temperature -dependence  to  the  exponential 
one,  would  give  no  noticeable  changes. 
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In  Figure  2,  the  total  magnetic  resistivity  for  the  different 
specimens  is  shown  as  a  function  of  temperature.  It  exhibits  two  im¬ 
portant  features:  First, the  magnetic  resistivity  shows,  similarly  to 
the  total  electrical  resistivity,  a  temperature  independent  contribu¬ 
tion  which  is  a  sensitive  function  of  the  material.  We  attribute  this 
residual  magnetic  resistivity  partly  to  domain  effects  and  partly  to 
magnetic  impurities.  Second,  the  remaining,  temperature-dependent 
part,  is  roughly  the  same  for  all  samples  and  is,  therefore,  attri¬ 
buted  to  the  scattering  by  spinwaves. 

o 

It  is  unfortunate  that  in  the  temperature  region  below  15  K  , 
which  is  the  most  easily  accessible  region  for  theory,  nothing  can  be 
said  about  the  electron- spinwave  scattering  for  the  following  two 
reasons.  Firstj  the  temperature  dependence  of  the  magnetic  resistiv¬ 
ity  is  approximately  within  the  experimental  accuracy,  and  second, 
the  magnetic  impurity  scattering  term  can  also  have  a  slight  temper¬ 
ature  dependence  and  can,  therefore,  no  longer  be  separated  from 
the  spin  wave  scattering  term. 

A  summary  of  the  results  is  given  in  Figure  3.  From  Figure  3a 

the  temperature  dependence  of  the  total  resistivity  is  seen  to  vary 
3  5  3 

between  a  T  *  law  and  a  T  law.  In  order  to  get  the  phonon  resistivi¬ 
ty,  we  have  to  subtract  the  magnetic  part,  which  is  plotted  in  Figure 
3b.  This  phonon  part  is  then  shown  in  Figure  3c.  Two  things  can  be 

seen  immediately.  The  temperature  dependence  is  not  at  all  a  Grtin- 

5  2  8  3  5 

eisen  T  law  but  varies  between  T  *  and  T  ’  .  Second,  in  our 

temperature  range,  the  phonon  part  is  considerably  larger  than  the 

magnetic  part,  contrary  to  what  is  obtained  by  subtracting  a  Griin- 
5 

eisen  T  law  from  the  total  resistivity,  A  temperature  dependence 

5 

different  from  a  T  law  can  be  expected  if  Umklapp-processes  are 
important  as  will  be  the  case  in  a  polyvalent  metal  such  as  Gd. 
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The  difference  in  magnitude  of  the  phonon  resistivity  in  various  sam¬ 
ples  is  due  to  different  textures  in  the  polycrystalline  samples,  which 
was  confirmed  by  X-ray  measurements. 

The  resistivity  arising  from  the  scattering  by  spinwaves  is 

plotted  in  Figure  3b.  With  the  exception  of  Gd  10,  the  spinwave  re- 

3.1 

sistivity  follows  roughly  a  T  *  law.  Gd  10  showed  an  extremely 
high  residual  magnetic  resistivity.  An  additional,  weak  temperature 
dependence  of  the  scattering  by  magnetic  impurities  and  domains 
could,  therefore,  lead  to  the  rather  large  temperature  dependent  part 
of  the  magnetic  scattering  at  low  temperatures.  The  difference  of  the 
spinwave  resistivity  for  various  amples  is  also  attributed  to  the  dif¬ 
ferent  texture  of  the  specimens. 

As  already  mentioned  above,  we  cannot  determine  the  magnetic 

resistivity  at  low  temperatures  accurately  enough  to  check  the  well- 
2 

known  Kasuya-T  law  valid  for  a  free  electron  gas.  Little  is  known 

2) 

of  the  band  structure  of  the  rare  earth  metals.  Recently,  Kasuya  ' 

suggested  a  model  which  leads  to  highly  anisotropic  energy  bands  and 

spinwave  spectra.  In  view  of  the  qualitative  character  of  this  model, 

it  is  too  early  to  predict  anything  about  the  magnetic  resistivity.  The 

only  quantitative  calculation  in  cur  temperature  range  was  made  by 
3) 

Goodings  /  which  is  shown  by  the  dotted  line.  He  considered  the  spin- 
wave  induced  s-d  scattering.  Considering  the  many  approximations 
he  uses,  the  agreement  or  disagreement  with  our  results  should  not 
be  taken  too  seriously. 

REFERENCES 

1.  M.  Kohler,  Ann.  Phys.  32,  211  (1938). 

2.  T.  Kasuya,  to  be  published  in  "Magnetism”  edited  by  G.  T,  Rado 
and  H.  Suhl. 

D.  A.  Goodings,  Phys.  Rev.  132,  542  (1963). 


3. 


THERMAL  CONDUCTIVITY  AND  ELECTRICAL  RESISTIVITY 
OF  SAMARIUM  AT  LOW  TEMPERATURES 
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ABSTRACT 

Thermal  conductivity  and  electrical  resistivity  of 
samarium  have  been  measured  from  5  to  200°K  and  from  5  to  300°K, 
respectively.  Two  anomalies,  one  at  12  +  1°K  and  another  at 
106  +  1°K,  are  observed,  resulting  from  an  order-order  magnetic 
transformation  and  an  antiferromagnetic-paramagnetic  transition, 
respectively.  The  Lorenz  function  is  found  to  be  larger  at  any 
temperature  than  that  expected  for  pure  electronic  thermal 
conductivity.  .This  implies  that  there  is  some  phonon  and 
possibly  also  some  magnon  thermal  conductivity  in  samarium  at 
low  temperatures.  The  magnetic  moment  disorder  electrical 
resistivity  of  samarium  is  determined  to  be  39.0  +0.5  pD  cm, 
in  fair  agreement  with  the  value  to  be  expected  from  theoretical 
considerations. 

Introduction 

According  to  studies  of  the  magnetic  susceptibility*, 
specific  heat^’^j  electrical  resistivity^  and  thermoelectric 
power  ,  polycrystalline  samarium  exhibits  two  (magnetic)  transitions 
at  low  temperatures.  One  occurs  at  about  12  +  1°K  and  the  other  at 
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106  +  1°K.  Neutron  diffraction  investigations  have  not  yet  been  done 

on  samarium.  Thus,  the  nature  of  these  transitions  is  not  clearly 

understood  at  the  present  time.  However,  on  the  basis  of  the  behavior 

of  the  above-mentioned  physical  properties,  it  appears  that  the 

transition  at  106°K  is  the  Nlel  point,  i.e.,  a  transition  from  an 

antiferromagnetic  state  to  the  paramagnetic  state.  The  low 

temperature  transition  is  quite  likely  a  magnetic  order-order 

8-i  i 

transition.  Recently  we  have  found  "  that  the  changes  in  the 

magnetic  state  in  gadolinium  dysprosium,  terbium,  and  erbium,  whose 

detailed  moment  arrangements  have  been  determined  by  the  neutron 

12-15 

diffraction  techniques  ,  exhibit  noticeable  influence  on  the 
heat  conductivity  of  these  metals.  Since  thermal  transport  studies 
have  not  yet  been  done  on  the  light  rare  earth  metals,  we  decided  to 
perform  such  measurements.  In  this  paper  we  present  the  results  of 
thermal  conductivity  and  electrical  resistivity  of  samarium  and 
discuss  their  significance. 

Experimental  Considerations 

The  sample  of  samarium,  whose  electrical  resistivity  and 
thermal  conductivity  were  measured  as  a  function  of  temperature,  was 
prepared  in  the  following  manner.  The  initial  stock  or  distilled 
samarium  was  obtained  from  Research  Chemicals.  According  to  the 
supplier  a  partial  analysis  showed  the  following  impurities  in  weight 
%i  Eu  -  0.0$,  Ca  -  0.02,  Mg  -  0.01,  Gd  -  0.01,  and  Si  -  0.005.  The 
metal  was  arc  melted  for  10  minutes  in  100  Torr  argon  atmosphere. 
After  the  melting,  a  rod  of  0.479  cm  diameter  and  about  6  cm  length 
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was  machined  from  the  ingot.  Swaging  of  samarium  ingots  resulted  in 
fracture,  and  it  was  therefore  impossible  to  use  this  method  to 
produce  a  good  polycrystalline  structure.  The  electrical  resistivity 
of  this  sample  at  4.17°K  was  found  to  be  6.73  pQ  cm.  The  electrical 
resistivity  and  thermal  conductivity  data  were  obtained  using  the 
equipment  described  elsewhere^’ 

Results  and  Discussion 

The  electrical  resistivity  of  samarium  from  liquid  helium 
temperatures  to  300°K  is  shown  in  Fig.  1.  One  clearly  can  observe 
two  anomalies  at  low  temperatures:  one  at  12°K  and  the  other  at  106°K. 
The  detailed  temperature  variation  of  the  electrical  resistivity  in 
the  neighborhood  of  these  temperatures  is  shown  in  Figs.  2  and  3, 
respectively.  Since  the  temperature  coefficient  of  the  electrical 
resistivity,  defined  by 

a  =  (l/p)(dp/dT),  (1) 

where  p  is  the  total  electrical  resistivity  at  a  temperature  T,  is  a 
particularly  sensitive  function  for  the  magnetic  transitions,  we  also 
have  plotted  this  quantity  in  Figs.  2  and  3.  From  these  plots  the 
above-given  transition  temperatures  were  determined. 

The  electrical  resistivity  of  samarium  changes  almost 
linearly  with  temperature  below  300°K.  At  about  106°K  the 
resistivity  decreases  much  more  rapidly  with  decreasing  temperatures 
than  above  this  temperature.  This  behavior  strongly  suggests  that 
above  106°K  the  magnetic  moments,  resulting  from  the  localized  4f 
electrons  in  samarium,  are  disordered  and  give  rise  to  considerable 
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moment  disorder  resistivity  .  We  estimate  this  from  the  experimental 

data  to  be  39.0  +  0.5  pfi  cm.  If  one  neglects  toe  transfer 'Of  energy 

between  the  conduction  electrons  and  the  magnetic  moment  system,  the 

20  21 

moment  disorder  electrical  resistivity  is  given  by  ’ 


p  j  =  3TT-Nry-  G  (g-l-)2  J(J+1), 
md  2  h  e2E_ 


(2) 


where  N  is  the  number  of  atoms  per  unit  volume,  m*  the  effective  mass 
of  a  conduction  electron,  e  the  electronic  charge,  Ep  the  Fermi  energy, 
g  the  Land4  g-factor,  J  the  quantum  number  associated  with  the  total 
angular  momentum  of  the  magnetic  moment,  and  h  =  h/2n,  h  being  the 
Planck  constant.  The  symbol  G  in  Equation  (2)  represents  a  coupling 
constant  with  the  dimensions  of  energy  times  volume.  Since  the 
quantity  in  front  of  (g-l)  J(J+l)  in  Equation  (2)  is  expected  to  be 
approximately  the  same  for  all  trivalent  rare  earth  metals,  one 
predicts  that  pm^  should  be  proportional  to  (g-l)  J(J+l).  Since 
this  factor  is  4.46  and  4.50  for  samarium  and  holmium,  respectively, 
resistivity  should  be  approximately  equal  for  these  two  rare  earth 

metals.  This,  indeed,  is  in  accordance  with  the  experimental  results: 

19 

for  holmium  p^  =  32.3  p£2  cm  which  should  be  compared  with  = 

39.0  pfi  cm  for  samarium. 

The  thermal  conductivity  of  samarium,  also  shown  in  Fig.  1, 
exhibits  two  maxima  at  low  temperatures:  one  at  about  9°K  and  another 
at  32oK.  The  first  maximum  and  the  following  minimum  at  13°K  results 
from  the  low  temperature  magnetic  order-order  transformation.  The 
observed  thermal  conductivity  behavior  is  qualitatively  consistent 


f*47^V 
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with  the  electrical  resistivity  behavior,  i.e.,  the  rapid  increase  in 
the  resistivity  above  10°K  causes  a  decrease  in  the  thermal  conduc¬ 
tivity.  The  second  maximum  is,  of  course,  of  the  conventional  type 
observed  in  all  pure  metals.  This  maximum,  qualitatively  can  be 
understood  if  one  recalls  that  at  low  temperatures  (T  «  0,  where  0 
is  the  Debye  temperature)  the  electronic  thermal  conductivity, 
limited  by  the  impurity  scattering,  is  proportional  to  T.  At  higher 

temperatures  (T  =»  0)  the  electron-phonon  scattering  causes  the 

_2 

electronic  thermal  conductivity  to  follow  approximately  T 
dependence.  Thus,  a  maximum  in  the  X  vs.  T  curve  results. 

The  thermal  conductivity  change  at  the  Nlel  point,  which 
according  to  the  electrical  resistivity  data  is  106°K,  is  not  abrupt 
but  gradual.  However,  the  general  increase  in  the  thermal  conduc¬ 
tivity,  when  the  magnetic  moment  becomes  disordered,  is  in  agreement 
with  the  results  on  the  other  rare  earth  metals^’ *  * ^  above  their 
N^el  temperatures. 

In  principle,  the  total  thermal  conductivity  of  a  rare 
earth  metal,  possessing  an  ordered  magnetic  structure,  should  be 
written  as 

X  =  Xe  +  X^  +  Xm  ,  (3) 

where  the  superscripts,  e,  £,  and  m  imply  that  the  heat  is  carried 
by  electrons,  phonons,  and  magnons,  respectively.  For  pure 
electronic  thermal  conductivity  at  low  temperatures  (T  «  0) ,  where 
the  scattering  processes  can  be  considered  elastic,  one  would  expect 
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that  the  Lorenz  function 

L  =  pX/T  (4) 

should  be  equal  to  the  Lorenz  number 

Lq  =  n2k2/3e2  =  2.445  x  10"8[V2°K“2].  (5) 

In  this  equation  k  is  the  Boltzmann  constant  and  e  the  electronic 
charge.  If  some  heat  is  also  transported  by  phonons  and  magnons, 
then  L  should  be  larger  than  Lq. 

Figure  4  shows  the  function  L  for  samarium  between  5  and 
200°K.  Although  the  Lorenz  function  was  not  determined  from 

simultaneous  measurements  of  p  and  X  (these’ properties  were  measured 

/ 

during  separate  runs  using  the  same  potential  contacts)  it  is  believed 

that  L  is  sufficiently  accurate  to  clearly  establish  that  the 

experimental  L  is  considerably  larger  than  Lq  at  any  temperature. 

This  implies  that  there  is  some  heat  conduction  by  phonons  and  also 

possibly  magnons.  This  aspect  of  the  quantity  L  for  samarium  is 

similar  to  the  behavior  found  in  dysprosium8,  gadolinium^,  terbium'1'8, 

and  erbium^.  Since  the  present  theoretical  understanding  of  the  role 

of  magnons  in  the  total  heat  conductivity  of  a  magnetically  ordered 

metal  is  quite  poor,  it  is  difficult  at  this  moment  to  separate  the 
^  rn 

terms  X  and  X  from  the  measured  values  of  X.  We  are  planning  in  the 
near  future  to  perform  thermal  conductivity  measurements  on  the  above- 
mentioned  rare  earth  metals  (and  also  other  ferromagnets)  in  the 
presence  of  strong  magnetic  fields  (up  to  60  KOe)  in  order  to  learn 
more  about  the  magnitude  and  temperature  dependence  of  Xm. 

Finally,  we  like  to  remark  that  recently  accurate  thermal 
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conductivity  measurements  on  rare  earth  metals  at  291°K  have  been  made 

22  -1 
by  Powell  and  Jolliffe  .  For  samarium  they  find  X  =  0.133  W  cm 

-1  2  -2 

°K  and  L  =  4.2  V  °K  .  These  values  appear  to  be  consistent  with 

our  measurements  when  extrapolated  from  200°K  to  291°K.  We  did  not 

perform  thermal  conductivity  studies  on  samarium  above  200°K  because 

of  serious  difficulties  experienced  with  radiation  connections.  For 

8—11 

this  reason  our  previous  data  on  other  rare  earth  metals  above 
liquid  nitrogen  temperatures  now  appear  to  be  higher  than  the  true 
values  as  can  be  concluded  by  the  work  of  Powell  and  Jolliffe. 
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ABSTRACT 

Ten  sets  of  electronic  band  structure  calculations  were  performed 
for  cerium  metal  by  using  the  augmented  plane  wave  method  for  five  con- 
figurations  of  the  type  6s  5<i  kf  and  the  lattice  constants  appropri¬ 
ate  to  the  FCC  gamma  and  FCC  alpha  phases .  It  is  shown  that  the 
location  of  the  f  band  is  strongly  dependent  on  the  atomic  configuration 
assumed  in  starting  the  calculation.  The  relation  of  these  results  to 
previous  atomic  calculations  is  discussed. 

It  is  believed  that  the  bands  should  not  be  occupied  to  a  common 
energy  level  according  to  a  simple  Fermi -Dirac  distribution,  but  rather 
that  the  4  electrons  outside  the  xenon  core  should  be  arranged  in  the  f 
as  *r  11  as  in  the  s  and  d  bands  so  as  to  be  reasonably  consistent  with 
the  configuration  assumed.  As  a  consequence  the  f-band  location  is 
dependent  on  the  degree  to  which  it  is  occupied.  The  two-electron  cal¬ 
culations  by  Waber,  Liberman  and  Cromer  of  the  total  energy  for  a  cerium 
atom  confined  to  a  Wigner-Seitz  sphere  lend  further  support  to  this 
assignment  hypothesis . 
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Introduction 

Cerium  metal,  because  of  its  several  allot ropic  forms,  offers  the 
interesting  opportunity  of  investigating  how-  the  electronic  band  struc¬ 
ture  of  a  metal  changes  vith  the  arrangement  and  spacing  cf  neighboring 
atoms.  The  band  structure  of  iron  in  the  body-centered  and  face- 
centered  cubic  forms  was  studied  by  Wood."^  Such  an  investigation  appar¬ 
ently  has  not  been  done  for  any  other  metal.  The  explanation  given  by 
2  ^ 

Pauling  and  by  Zachariasen^  for  the  transformation  of  cerium  from  the 

FCC  gamma  form  to  the  "collapsed"  FCC  alpha  .form  is  that  a  4f  electron 

is  promoted  into  the  conduction  band  when  the  alpha  phase  is  formed. 

To  date  the  only  study  of  the  effect  of  configuration  on  band  structure 

4 

was  done  by  Mattheiss  for  vanadium.  Even  though  only  a  portion  of  the 
bands  was  calculated,  sensitivity  to  assumed  configuration  was  indicated. 
For  tir'e  reasons  cerium  is  an  interesting  metal  to  study  in  seme  detail. 
The  baau  structures  of  the  double  HOP  beta  and  the  BCC  delta  phase  of 
cerium  will  not  be  discussed. 

The  only  rare-earth  metal  for  which  a  band  calculation  has  been 

5 

done  to  date  is  gadolinium,  studied  by  Dimmock  and  Freeman.  They 

7  12 

assumed  the  configuration  to  be  4f ’5d  6s  .  They  found  that  the  narrow 
f  bands  were  approximately  0.8  Rydberg  below  the  s  and  d  bands.  Rather 
than  permitting  the  three  conduction  electrons  to  occupy  this  f  band, 
they  argued  that  to  regard  gadolinium  as  being  in  a  configuration  f^ 
would  substantially  raise  the  total  energy.  Thus  from  their  experience 
with  Hartree-Fock  calculations  (which  are  of  a  two-electron  nature)  they 
chose  a  configuration  before  solving  the  one-electron  band  structure 
problem.  Apparently,  however,  no  numbers  were  presented  to  justify 
their  plausible  approach. 

6-9 

Raynor  and  his  colleagues  have  attempted  to  elucidate  the  change 
in  the  metallic  valence  cf  cerium  dissolved  in  different  metals.  Assum¬ 
ing  that  the  valence  of  a  metal  can  be  directly  related  to  its  metallic 
diameter  they  deduce, from  the  apparent  atonic  diameter,  that  cerium 
acquires  (or  approaches)  the  valence  of  the  solvent  metal.  The  opera¬ 
tive  mechanism  they  postulate  is  the  variation  of  strain  energy  with 


77 


the  size  of  the  hole  into  which  the  solute  is  placed. 

No  proposed  electronic  configuration  for  cerium  metal  has  been  uni¬ 
versally  accepted.  Gschneidner  and  Smoluchowski  concluded  on  the 
basis  of  a  variety  of  experimental  data  that  one  4f  electron  is  in  a 
localized  atomic  state  in  gamma  cerium  and  approximately  O.k  of  a  4f 
electron  in  the  alpha  phase.  In  the  present  study,  no  a  priori  division 
of  electrons  was  made.  Certain  conclusions  can  nevertheless  be  made. 

Previous  Studies  of  Atomic  Configuration 
11 

Larson  and  Waber  presented  the  results  of  one-electron  Hartree 
calculations  for  the  free  cerium  atom  at  the  Third  Rare  Earth  Conference. 
The  binding  energy  of  the  4f  electrons  in  the  general  atomic  configura- 
tion  6s  5<i  ^f  *  was  shown  to  increase  with  the  occupancy  of  the  5d 
orbitals ,  From  the  total  energy  of  the  free  atom  in  this  Hartree  approx¬ 
imation  they  concluded  that  the  probable  configuration  of  the  free  atom 
was  6s25d1'^f0'^;  the  total  energy  was  approximately  -16,79*4-  Rydbergs. 

12 

As  a  more  reasonable  approximation,  Waber,  Liberman  and  Crcmer 

studied  for  the  Fourth  Rare  Earth  Conference  the  behavior  of  a  cerium 

atom  with  the  same  general,  configuration  confined  to  a  Wigner-Seitz 

sphere  of  radius  R^.  In  place  of  the  non-relativistic  self -consistent 

field  (SCF)  Hartree  functions,  they  used  the  relativistic  SCF  Dirac- 

13 

Slater  wave  functions,  x  which  include  exchange.  The  strong  dependence 
of  the  one -electron  energy  eigenvalues  on  the  occupancy  of  the  5&  orbi¬ 
tals  is  shown  in  Fig.  1.  Curves  are  shown  for  two  radii,  namely  the 
Wigner-Seitz  radii  appropriate  to  the  alpha  (dotted)  and  gamma  (full) 
phases . 

This  graph  leaves  the  impression  that  the  most  stable  configuration 
2  2 

is  6s  5<i  )  since  each  of  the  eigenvalues  is  more  negative  for  x  =2.0. 
When  two-electron  effects  are  included,  this  conclusion  is  modified.  The 

ik 

total  energy  E^  was  computed  by  the  method  of  Snow,  Canfield  and  Waber. 
Results  for  different  radii  R^  and  different  configurations  are  shown 
in  Fig.  2.  One  concludes  that  the  most  stable  configuration  depends  on 
the  size  of  the  Wigner-Seitz  sphere,  i.e.  on  the  pressure. 
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It  will  be  shown  below  that  in  the  band  structure  calculations  the 
location  of  the  bands  depends  on  the  assumed  electronic  configuration. 

The  f-like  bands  depend  strongly  on  x,  the  occupancy  of  the  d  bands . 

Method  of  Calculation 

The  augmented  plane  wave  method  was  first  proposed  in  1937  by 
15 

Slater.  The  principal  computational  details  of  this  method  have  been 
given  by  Wood."'"  In  the  present  work,  the  Bloch  function  was  expanded  in 
terms  of  spherical  harmonics  with  angular  momenta  Z  <;  13.  In  general, 
the  maximum  value  of  the  wave  vectors  retained  in  the  set  was  6n/aQ. 

With  this  set,  convergence  was  indicated  to  be  better  then  0.003 
Rydberg* 

'\Jo 

Em  and  Switendick  have  discussed  how  to  obtain  the  charge  den¬ 
sity  associated  with  a  particular  angular  momentum  from  the  symmetrized 
AFW  crystal  eigenfunctions. 

17 

The  crystal  potential  was  constructed  by  the  method  of  Mattheis.  1 

Namely,  atomic  overlap  contributions  were  computed  by  Lowdin's  alpha 
l8 

expansion  method.  Contributions  to  the  electronic  charge  densities 
were  added  out  to  the  set  of  fourth  nearest  neighbors,  in  the  present 
study,  and  spherically  averaged.  Dirac -Slater  wave  functions  were  used 
for  each  configuration.  Slater’s  method"^  was  employed  to  incor¬ 
porate  exchange  by  using  the  summed  and  spherically  averaged  charge 
density  P- 

Since  the  AFW  program  available  to  us  used  the  "telescoping"  linear 

20 

presentation  mesh  of  Herman  and  Skillman,  values  of  the  potential  were 

interpolated  at  appropriate  radii  from  the  logarithmic  radial  mesh  that 

13 

Liberman,  Waber  and  Cromer  used  in  calculating  the  Dirac -Slater  wave 
functions. 

The  constant  or  muffin  tin  potential  in  the  region  between  the  AFW 
sphere  and  the  Wigner-Seitz  sphere  was  obtained  by  spherically  averaging 
the  summed  potentials  and  the  exchange  contribution.  It  was  found  that 
this  potential  increased  linearly  with  x,  the  occupancy,  in  both  the 
gamma  and  alpha  size  spheres.  The  numbers  are  tabulated  in  Figs.  3  and  4. 
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Results 

The  band  structure  curves  have  been  obtained  for  two  sizes  of  the 

lattice  constant  and  for  five  electronic  configurations  in  each  size. 

The  complicated  energy  dependence  E(k)  was  obtained  at  256  points  in  the 

Brillouin  zone,  namely  at  points  ^  >  ?)  0 f  the  reciprocal  lattice. 

d  d  d 

Because  of  the  narrowness  of  the  f  bands,  it  was  necessary  to  evaluate 
the  secular  determinant  at  energy  intervals  as  small  as  8  x  10 Rydberg 
or  9  cm  The  results  reported  here  may  be  regarded  as  preliminary 
since  it  was  difficult  to  obtain  all  13  eigenvalues  (for  s,  d  and  f-like 
bands)  associated  with  certain  low  symmetry  or  general  wave  vectors. 

The  over-all  behavior  of  the  bands  as  a  function  of  configuration 
is  shown  in  Figs.  3  and  4.  For  convenience  and  clarity,  only  the  direc¬ 
tion  T  to  X  (along  A)  is  shown  for  each  graph.  The  AFW  method  partially 
incorporates  the  effect  of  the  crystal  field  in  lifting  the  degeneracy 
of  states  having  a  given  angular  momentum.  For  example,  the  seven  f- 
like  atonic  states  (either  spin)  are  separated  into  a  non-degenerate 
Fg^  and  the  two  triply  degenerate  ^  and  states.  Similarly,  the 
5d  atomic  states  are  split  into  the  triply  degenerate  (tg  )  and 

doubly  degenerate  ^(e^)  states.  Comparing  Fig.  3  for  the  gamma- phase 
size  with  Fig.  4  for  the  collapsed  or  alpha  phase,  one  sees  that  the 
crystal  field  splitting  is  larger  in  the  latter  case.  In  addition,  one 
observes  that  the  band  width  is  larger,  presumably  because  of  increased 
overlap. 

The  pattern  of  the  s  and  d  bands,  such  as  is  typical  of  a  3d  transi¬ 
tion  metal,  is  clearly  discernible  at  the  far  right-hand  side  of  these 
figures  where  it  is  assumed  that  no  4f  states  are  occupied.  This  premise 
is  unjustified  because  the  calculations  yield  low- lying  f  states.  At 
the  far  left-hand  side,  the  f  bands,  which  are  shown  as  a  cross-hatched 
region,  interact  with  the  s  and  d  bands  and  the  result  is  a  more 


The  notation  used  is  substantially  the  same  as  that  of  Bouckaert, 
Smoluchowski  and  Wigner^  except  that  the  primed  states  are  designated 
by  asterisks  to  facilitate  reading  of  the  graphs. 
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2  2 

complicated  picture.  Although  the  configuration  6s  4f  is  assumed,  the 
f  states  are  too  high  and  would  not  be  occupied;  that  is,  the  lower 
lying  states  which  will  hold  four  electrons  have  predominantly  s  eri  d 
character.  For  intermediate  configurations  such  as  one, 

on  the  basis  of  the  usual  one-electron  band  picture,  would  put  the  four 
available  electrons  into  the  low-lying  f  states,  despite  the  assumed 
configuration  with  a  large  d-character.  Thus,  one  is  led  to  conclude 
that  an  intermediate  value  of  x,  the  d- orbital  occupancy,  is  more 
reasonable . 


The  Schrodinger  equation,  as  solved  in  any  band  structure  calcula¬ 
tion  made  to  date,  involves  an  average  one-electron  potential.  That  is, 
all  solutions  are  obtained  in  the  field  of  a  single  assumed  charge  dis¬ 
tribution.  One  calculates  E(k)  curves  which  show  crystal  eigenstates 
of  the  assumed  potential.  These  states  are  occupied  by  the  valence 
electrons  in  accordance  with  Fermi -Dirac  statistics,  and  occupancies  are 
usually  found  which  are  reasonably  consistent  with  the  assumed  configura¬ 
tion.  In  contrast,  if  one  finds,  as  we  do  in  cerium,  a  derived  occupancy 
which  is  in  conflict  with  the  assumed  one,  other  methods  of  assigning  the 
occupancy  should  be  studied.  The  simple  Fermi -Dirac  distribution  may 
have  to  be  modified.  The  reason  is  that,  in  the  latter  model.,  the 
energy  states  are  assumed  to  be  independent  of  the  occupancy  whereas  our 
results  clearly  show  the  opposite. 


Wien  one  calculates  the  E(k)  curves  for  states  lying  above  the 
Fermi  energy,  what  he  does  in  analogy  to  the  atonic  SCF  calculations  is 
to  compute  the  eigenvalue  for  an  unoccupied  orbital.  This  often  is  a 
reasonable  approximation  to  the  correct  eigenvalue.  However,  in  a  more 
correct  treatment,  if  he  were  to  take  an  electron  frcm  a  normal  valence 
state  and  put  it  into  an  excited  orbital  such  as5g  and  then  complete  an 
SCF  calculation,  he  would  find  a  different  potential  due  to  the  new 
charge  distribution.  Consequently  all  of  the  one-electron  eigenvalues 
would  be  changed.  The  Is  electrons  would  be  affected  only  slightly  and 
their  charge  distribution  would  be  virtually  unaltered.  However,  for  the 
valence  electrons,  both  their  charge  distribution  and  their  eigenvalues 


would  be  strongly  affected.  The  rearrangement  energy  or  shifts  in 
eigenvalues  corresponding  to  occupancy  of  higher  states ,  or  states  not 
initially  assumed,  are  not  usually  considered  in  band  structure  calcula¬ 
tions.  The  present  study,  as  do  the  atomic  calculations,  shows  that  the 
rearrangement  energies  due  to  reassignment  of  the  occupancy  need  not  be 
small.  For  cerium,  rearrangement  energies  play  an  important  role.  A 

similar  situation,  although  not  as  drastic,  has  been  found  in  the  study 
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of  the  band  structure  of  nickel  by  Snow,  Svitendick  and  Waber. 

When  empty,  the  f  states  are  very  well  bound  as  Fig.  1  shows.  How¬ 
ever,  when  they  are  filled  to  a  reasonable  extent,  the  electrostatic 
repulsion  between  orbitals  increases  sharply,  and  their  eigenvalues  rise 
above  those  of  possible  competing  s  and  d  states.  Thus,  it  is  suggested 
that  a  balance  occurs,  which  results  in  an  intermediate  number  of  occu¬ 
pied  states.  These  general  remarks  apply  to  band  as  well  as  to  atcmic 
calculations. 

Generally  the  shape  of  the  E(k)  craves  is  relatively  insensitive  to 
the  assumed  electronic  configuration,  but  this  need  not  always  be  true. 

In  the  case  of  cerium.,  as  the  f  bands  move  through  the  region  of  the  s 
and  d  bands  considerable  modification  occurs  because  of  interaction 
between  bands. 

23 

Kittel  in  his  appendix  N  shows  that  the  partial  molar  free  ener¬ 
gies  (or  chemical  potentials  p,)  of  transferring  electrons  from  phase  i 
to  phase  j  must  be  equal.  In  equilibrium  the  Fermi  level  must  be  the 
same  in  the  two  phases  in  contact.  This  argument  can  easily  be  extended 
to  two  competing  bands.  The  partial  molar  free  energy  or  the  chemical 
potential  for  transferring  an  electron  from  the  s  and  d  bands  to  the  f 
bands,  or  vice  versa,  also  involves  the  energy  shift  due  to  the  rearrange¬ 
ment  of  occupancies .  Thus  the  energy  which  would  be  gained  by  the  trans¬ 
fer  of  an  electron  into  a  partially  filled  f  band  is  partially  offset  by 
the  increased  electrostatic  interactions  and  the  relocation  of  the  f 
band.  Such  interactions  are  not  embraced  in  the  usual  framework  of 
Fermi -Dirac  statistics.  Thus  we  believe  that  each  band  does  not  have  to 
be  occupied  to  the  same  level,  but  rather  that  partial  molar  free 
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energies  must  be  equal.  Because  of  the  effect  of  occupancy  in  the  f 
bands  on  their  location,  we  believe  that  the  most  appropriate  band  cal¬ 
culation  to  be  used  in  estimating  physical  properties  is  one  which  is 
most  nearly  consistent  with  the  assumed  configuration  or  occupancy.  In 
the  present  calculations,  intermediate  values  of  x  between  1.0  and  1.5 
in  the  general  configuration  6s  5d  4f  appear  to  meet  this  criterion. 

Certain  trends  are  clearly  indicated  by  the  present  study.  One  of 
these  is  the  effect  of  pressure.  If  we  look  at  the  6s  5d  4f  configura¬ 
tion  for  both  the  gamma  and  alpha  phases,  we  see  that,  for  alpha  cerium, 
the  f  bands  are  contiguous  to  the  s  and  d  bands,  whereas  in  gamma  cerium 
about  0.1  Rydberg  separates  the  f  band  from  the  s  and  d  bands.  Despite 
being  labeled  as  f  states,  there  is  considerable  mixing  and  a  consider- 
able  amount  of  d-character  when  they  are  near  the  d  band.  Thus  one 
would  conclude  that  a  consistent  assignment  of  occupancy  would  include 
more  d-like  electrons  in  the  collapsed  alpha  phase  than  in  the  gamma. 

Hie  present  results  indicating  an  inteimediate  value  of  x,  the  d 
orbital  occupancy,  and  a  larger  value  of  x  in  the  alpha  phase  are  in 
reasonable  agreement  with  the  two-electron  cellular  calculations  of  Fig. 
2.  Further,  the  results  are  not  inconsistent  with  the  assignment  which 

Gschneidner  and  Smoluchowski^0  made,  namely  for  gamma 

j  c.  2Cj1«5^'),.£>0  .^6  „  it, 
and  6s  5d  4f  for  alpha  cerium. 

One  of  the  inferences  of  this  investigation  is  that  a  satisfactory 
and  self-consistent  band  structure  of  a  general  metal  cannot  be  achieved 
until  the  two-electron  interactions  are  incorporated  into  the  analysis. 

As  a  first  step,  it  appears  to  be  reasonable  to  calculate  the  electro¬ 
static  <  ij|- — — ~rH  >  ssfi  the  exchange  <  ijj- — - j  ji  >  contribu- 

1  2  i  ‘  r2 

tions.  It  is  an  interesting  point  that  use  of  the  one-electron  API'/ 

method  yields  results  similar  to  those  found  by  the  two-electron  total 
energy  calculation,  namely  an  intermediate  value  of  x  is  a  more  consistent 


On  the  basis  of  our  charge  density  calculations  we  find  less  than  two 
s-like  electrons  in  the  conduction  band.  This  result  is  also  usually 
found  for  the  transition  metals. 


picture  of  occupancy  of  cerium  in  the  solid. 

To  sum  up,  a  large  but  reasonable  sensitivity  of  the  location  of 
S;  d  and  f  bands  to  electronic  configuration  assumed  for  cerium  atoms 
has  been  demonstrated.  Because  of  new  questions  raised  by  this  result, 
it  is  premature  to  assign  a  definite  configuration  to  metallic  cerium 
from  these  ab  initio  calculations. 
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TOTAL  ENERGY  +17700  R  Y 


VARIATIONS  OF  TOTAL  ENERGY  WITH 
CONFIGURATION  AND  SPHERE  SIZE 


X  NUMBER  OF  5d  ELECTRONS 


Fig.  2.  Effect  of  assumed  configuration  and  sphere  size  on  the  two- 
electron  total  energy  of  the  confined  cerium  atom. 


EFFECT  OF  ATOMC  CONFIGURATION  ON  THE  BAND  STRUCTURE  OF  GAMMA  CO*UM 


5.  Effect  of  the  assumed,  configuration  on  the  band  structure  of 
gamma  cerium. 


EFFECT  OF  ATOMK  CONFIGURATION  ON  THE  BANO  STRUCTURE  OF  ALPHA  CUWJW 


4.  Effect  of  the  assumed  configuration  on  the  band  structure  of 
alpha  cerium. 
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THEORETICAL  RESULTS  FOR  POSITRON  ANNIHILATION 
IN  RARE  EARTH  METALS* 

T.  L.  Loucks 

Institute  for  Atomic  Research  and  Department  of  Physics 
Iowa  State  University,  Ames,  Iowa 

ABSTRACT 

The  angular  correlation  of  the  photons  from  two-quantum  positron 
annihilation  has  recently  been  measured  in  single  crystals  of  holmium 
at  300°K  and  70°K.  ^  These  results  exhibit  some  very  interesting  struc¬ 
ture  as  compared  to  the  usual  free-electron  parabola,  and  they  show 
clearly  that  the  Fermi  surface  of  holmium  is  highly  anisotropic.  An 
attempt  is  being  made  to  reproduce  these  results  theoretically;  the 
method  of  calculation  will  be  discussed  in  this  paper. 

The  model  is  based  on  the  muffin- tin  potential  and  Augmented 
Plane  Wave  method  (APW)  of  Slater.  ^  The  potential  for  the  electrons 
(including  free-electron  exchange  based  on  the  local  charge  density  in 
the  crystal)  is  constructed  from  self-consistent- field  atomic  calculations. 
A  similar  potential  (with  opposite  sign)  is  constructed  for  the  positron 
with  exchange  excluded.  The  wave  function  for  the  positron  is  taken  to 
be  constant  outside  the  APW.  spheres  and  is  determined  numerically  in¬ 
side  the  spheres  using  the  positron  potential  and  requiring  zero  slope 
at  the  sphere  radius.  The  electronic  wave  function  is  taken  to  be  a 
linear  combination  of  APW's.  The  expansion  coefficients  are  deter¬ 
mined  by  finding  the  crystal  eigenvalues  from  the  secular  equation  and 
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solving  the  coupled  equations  which  result  from  the  variational  expres¬ 
sion  for  the  energy. 

The  probability  that  a  photon  pair  will  be  produced  with  z- compo¬ 
nent  of  momentum  Pz  is  proportional  to 

00  00 

N<Pz>  =  I  1  dPxdPy  P<Px-PyP.> 

00  _co  ' 

where 

p(P)  =  s  I  J  *k(r)  ^+(r)riP'~dr| 

k  — 

An  expression  for  N(Pz)  has  been  developed  using  the  above  muffin-tin 
model  for  the  electron  and  positron  wave  functions.  Since  it  is  based 
on  the  APW  method  it  will  have  application  not  only  to  free- electron 
metals  but  also  to  transition  and  rare-earth  metals.  The  results  of  this 
method  should  be  very  useful  as  an  aid  in  the  interpretation  of  positron 
annihilation  experiments  since  these  depend  in  a  very  complicated  way 
on  both  the  topology  of  the  Fermi  surface  and  the  wave  functions  of  the 
annihilating  particles. 


^Work  performed  in  the  Ames  Laboratory  of  the  U.  S.  A.  E.  C. 

1.  A.  R.  Mackintosh  and  R.  W.  Williams  (to  be  published). 

2.  J.  C.  Slater,  Phys.  Rev.  51,  846  (1937). 
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HIGH  PRESSURE  POLYMORPHISM  AND  MAGNETIC  ORDERING  IN 
SOME  Sm-TYPE  RARE  EARTH  ALLOYS  AND  SAMARIUM 

A.  Jayaraman,  R.  C.  Sherwood.,  and  E.  Corenzwit 

Bell  Telephone  Laboratories,  Incorporated 
Murray  Hill,  New  Jersey 

ABSTRACT 

The  rare  earth  alloys  LaY,  PrY,  NdY,  NdTm  and  LaGd^ 
all  of  which  crystallize  in  the  Sm-type  structure  undergo 
a  pressure-induced  phase  change  to  the  double  hexagonal 
close  packed  structure  (La-type  )j  the  high  pressure  phase  is 
metastably  retained  at  atmospheric  pressure.  However, 
under  similar  pressure  and  temperature  treatment  CeY  decom¬ 
poses  to  the  elements.  Magnetic  data  are  presented  for  the 
normal  and  high  pressure  phases  of  NdTm,  LaGd^  and  Sm;  in 
the  latter  case  the  data  are  from  earlier  measurements. 

The  significance  of  the  breakdown  of  CeY  and  the  anomalies 
in  the  low  temperature  magnetization  data  of  LaGd^,  NdTm 
and  Sm  are  discussed. 

INTRODUCTION 

A  pressure-induced  phase  transformation  to  the  double 
hexagonal  close  packed  structure  has  been  reported  in  Sm- 
metal.1)  The  high  pressure  phase  (dhcp)  is  retained  meta¬ 
stably  at  atmospheric  pressure  and  undergoes  antiferromag¬ 
netic  ordering  at  27°K,  about  12°  higher  than  the  Neel 
temperature  observed  for  the  normal  form  of  Sm.  Many  rare 
earth  binary  alloys  involving  a  lighter  rare  earth  element 
than  Sm,  and  a  heavier  rare  earth  element  or  Y,  crystallize 

p  o  li 

in  the  Sm-type  structure,  ,Di  near  the  equiatomic 
proportion.  In  some  systems  as  for  instance  La-Gd,  the 


i 


Sm-type  structure  is  stable  near  the  composition  LaGd^. 

It  occurred  to  us  that  these  alloys  might  undergo  a  pres¬ 
sure-induced  structural  transformation,  similar  to  that 
found  in  the  Sm  metal.  We  have  investigated  LaY,  CeY, 

PrY,  NdY,  NdTm  and  LaGd^  in  the  present  study.  The  mag¬ 
netic  behavior  of  NdTm  and  LaGd^  was  studied  both  in  the 
normal  as  well  as  in  the  high  pressure  form.  The  results 
will  be  presented  and  discussed  in  this  paper. 

EXPERIMENTAL 

The  alloys  were  prepared  by  arc-melting  in  an  argon 
atmosphere  appropriate  amounts  of  the  metals.  They  were 
then  heat  treated  at  about  500°C  for  several  days,  in 
sealed  evacuated  quartz  glass  tubes.  X-ray  powder  patterns 
revealed  that  the  annealed  samples  were  single  phase  mate¬ 
rials  and  had  the  Sm-type  structure.  High  pressure  and 
temperature  treatment  was  carried  out  in  a  piston  cylinder 
apparatus.  Samples  were  normally  subjected  to  40  kbar  and 
500 °C  for  about  5  hours.  After  this  treatment  temperature 
was  reduced  to  ambient  and  pressure  released.  The  X-ray 
patterns  of  the  specimens  were  then  recorded.  Magnetic 
susceptibility  measurements  were  carried  out  with  a  pendu¬ 
lum  magnetometer^  down  to  1.4°K. 

RESULTS 

All  the  alloys  except  CeY  transformed  to  dhcp  struc¬ 
ture  (La-type)  and  the  high  pressure  phase  is  retained  on 
release  of  pressure.  Under  similar  pressure  and  tempera¬ 
ture  treatment,  CeY  broke  down  to  the  elements.  In  Table  I 
the  pertinent  crystallographic  data  are  presented.  Al¬ 
though  no  attempt  was  made  to  locate  the  transition 
pressures  precisely,  it  would  seem  from  the  results  that 
the  transition  pressure  increases  in  the  sequence 
LaY  <  PrY  <  NdY,  With  NdY  conversion  to  dhcp  was  only 


Table  I.  Crystallographic  Data 
Sm-Type  Phase  La-Type  Phase 

a(A)  c(A)  a(A)  c(A) 

LaGd3  3.65±.02  26.33±.01  3.63+.02  ll.79i.01 

LaY  3.701.02  26. 651. 01  3.691,02  ll.85i.Ol 

PrY  3,651.02  26.251.01  3.621.02  11. 681.01 

NdTm  3.601.02  26. 001. 01  3.601.02  11. 501.01 

partial  under  the  experimental  conditions  stated  above. 
Preliminary  studies  on  NdHo  and  NdGd  which  have  also  the 
Sm-type  structure  indicate  that  a  transformation  to  the 
dhcp  structure  (La-type)  takes  place  in  these. 

The  results  of  magnetization  measurements  are  pre¬ 
sented  in  Pigs.  1  (Sm-form)  and  2  (dhcp-form)  for  LaGd^ 
and  Pigs.  3  (Sm-type)  and  4  (dhcp- type)  for  NdTm.  In  the 
case  of  LaGd^  a  kink  in  the  curve  occurs  at  about  120°K 
for  the  Sm-form  (See  Pig.  l).  The  curve  for  NdTm  shows  a 
peak  for  the  Sm-form  at  about  32 °K  (See  Fig.  3).  In 
either  case,  the  curves  for  the  dhcp  form  appear  to  be 
smooth.  Pig.  5  is  the  data  for  Sm  taken  at  15300  oersteds 
and  represents  our  earlier  measurements.^  Curve  1  is  the 
data  for  the  dhcp  phase  and  2  for  the  normal  form  of  Sm. 

The  latter  has  a  kink  in  the  curve  at  about  110 °K  and  the 
dhcp  form  shows  a  very  much  enhanced  anomaly  at  about 
100°K.  A  sharp  peak  in  magnetization  occurs  at  l4.8°K  for 
the  normal  phase  of  Sm  and  at  28°K  for  the  dhcp  phase. 

DISCUSSION 

Crystal  structure  and  size  relationship 

In  the  rare  earth  metals  pressure  induced  phase  trans¬ 
formations  have  been  shown  to  occur  in  the  sequence  hep  -» 
Sm-type  -»  dhcp  -»  fee.'1')  The  results  of  the  present  study 
indicate  that  the  inura-rare  earth  alloys  also  fit  into 
this  sequence;  the  Sm-type  invariably  transforming  to  the 
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Magnetic  data  for  Sm-form  of  NdTm. 


Magnetic  data  for  dhcp  NdTm. 
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Pig.  5  Magnetic  data  for  Samarium. 
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dhcp  structure  (La-type)  under  pressure. 

The  behavior  of  CeY  can  be  rationalized  as  follows. 
Gschneidner  et  al.^  have  investigated  the  effect  of  the 
addition  of  Y  on  the  y  to  a-Ce  transformation  and  find 
that  the  transition  pressure  rises  with  the  addition  of  Y. 
In  the  y  to  a-Ce  transition  the  metallic  radius  of  Ce 
changes  from  1.824  A  to  1.71  A.10)  It  is  suggested  that 
the  Ce  atoms  in  CeY  undergo  this  reduction  in  radius 
under  pressure,  resulting  in  a  drastic  change  in  the  Ce  to 
Y  radius  ratio.  It  must  be  that  this  change  in  radius 
ratio  is  incompatible  with  the  stability  of  the  Sm-type 
structure  and  causes  the  breakdown  of  the  Sm-type  CeY  to 
the  elements .  In  Table  II  are  tabulated  the  radius  ratio 

Table  II.  Radius  Ratio  For  the  R.E.  Alloys 


La 

Ce3+ 

Ce4+ 

Pr 

Nd 

Y 

1.06 

1.026 

0.956 

1.032 

1.028 

Tm 

- 

- 

- 

1.034 

Ho 

1.053 

1.02 

- 

- 

1.022 

Gd 

i.o4i 

1.008 

— 

1.01 

(metallic  radius)  of  the  lighter  rare  earth  to  the  heavier 

R«E.  element  or  Y  for  several  alloys  crystallizing  in  the 

Sm-type  structure.  For  Gd  which  transforms  under  high 

pressure  and  temperature  (~35  kbar  and  350°C)  to  the  Sm- 

type  structure'1'1)  and  for  Sm,  this  ratio  is  unity.  In  the 
44- 

case  of  Ce  the  radius  ratio  would  be  substantially  lower 
than  unity  and  is  the  only  exception  (See  Table  II).  Thus 
it  appears  that  the  stability  of  the  Sm-type  structure  in 
the  intra-rare  earth  alloys  is  controlled  by  the  size  of 
the  constituent  atoms  and  apparently  for  a  radius  ratio 
less  than  unity,  neither  the  Sm-type  structure  nor  the 
dhcp  structure  is  stable. 
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Magnetic  behavior 

The  magnetic  behavior  of  Gd-La  and  Gd-Y  alloys  has 

12  \ 

been  studied  by  Thoburn,  Legvold  and  Spedding.  '  Accord¬ 
ing  to  these  authors  LaGd^  has  a  poorly  defined  antiferro¬ 
magnetic  maximum  at  about  130°K.  In  our  measurements  on 
LaGd^  (Sm-type),  an  anomaly  in  the  magnetization  curve  is 
evident  at  about  120°K.  The  Sm  form  undergoes  probably 
some  kind  of  ordering  at  this  temperature.  In  the  magnet¬ 
ization  curve  for  the  dhcp  form  any  such  anomaly  is  not 
evident.  The  effective  moment  per  Gd  calculated  from  the 
susceptibility  data  in  the  paramagnetic  region  yields  8.45 
and  8.8  Bohr  magnetons  for  the  Sm  and  dhcp  forms  of  LaGd^ 
respectively.  At  4.2°K  the  magnetization  increases  with 
the  field  and  at  80  kilo  gauss,  Gd  is  contributing  1.82 
and  1.62  Bohr  magnetons  to  the  moment  in  the  Sm-form  and 
dhcp  form  respectively  and  the  magnetization  curves  show 
very  definitely  some  remanence  in  both  forms.  It  is  there¬ 
fore  believed  that  at  very  low  temperatures  (4.2°K)  both 
forms  are  in  all  probability  ferrimagnetic .  Another 
possibility  is  that  a  very  large  crystal  anisotropy  pre¬ 
vents  the  attainment  of  a  higher  moment  and  in  both  forms, 
LaGd^  is  ferromagnetic.  Lacking  single  crystal  data  we 
are  unable  to  be  more  definite  on  this .  In  the  Gd-Y 
alloys12^  an  alloy  of  composition  95 %  Gd  and  25$  Y  is 
ferromagnetic  and  Gd  contributes  its  full  moment  7*63  M-g 
at  a  field  of  17*5  k-oersteds.  However  the  Gd-Y  alloys 
are  hep  in  the  entire  compositional  range  and  apparently 
have  low  crystal  anisotropy. 

13 ) 

Davis  and  Bozorth  '  studied  Tm  metal  and  reported  a 
maximum  in  the  susceptibility  curve  indicative  of  anti¬ 
ferromagnetic  ordering  at  about  6o°K.  According  to  them 

Tm  should  be  ferromagnetic  below  22 °K  as  it  shows  a  hys- 

14) 

teresis  loop.  Jennings,  Hill  and  Spedding  '  have  found 
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&  specific  heat  anomaly  at  55 °K  in  Tm  metal,  presumably 
connected  with  magnetic  ordering.  Koehler,  Cable,  Wollan 
and  Wilkinson1^)  studied  by  neutron  diffraction  the  mag¬ 
netic  structure  of  Tm  down  to  very  low  temperatures  and 
have  reported  evidence  for  an  oscillatory  type  antiferro¬ 
magnetic  configuration  (spiral  or  helical  type)  below 
55°K.  The  structure  adopted  by  Tm  at  very  low  tempera¬ 
tures  is  a  type  of  antiphase  domain-type  structure  in 
which  there  is  a  net  moment  of  1.0  Bohr  magneton,  due  to 
incomplete  cancellation.  For  an  alloy  of  89.4$  Tm-10.6^ 
La"^  which  is  hexagonal  close  packed,  Koehler  et  al. 
found  that  the  Neel  temperature  lowered  to  36 °K  with  the 
Curie  point  at  31°K.  NdTm  in  the  Sm-type  structure  would 
appear  to  be  similar  to  Tm  in  magnetic  behavior  in  many 
respects.  From  the  anomaly  in  the  magnetization  curve 
(See  Fig.  3)  it  is  deduced  that  NdTm  (Sm-type)  orders 
antiferromagnetically  at  about  30°K.  It  exhibits  a  hys¬ 
teresis  loop  at  very  low  temperatures  (our  data  was  taken 
at  1.3°K)  and  is  presumably  therefore  ferro  or  ferrimag- 
netic.  The  Curie  temperature  would  seem  to  be  near  about 
10°K.  There  does  not  appear  to  be  any  susceptibility 


anomaly  indicative  of  ordering  in  the  dhcp  form  of  NdTm 
near  30°K  and  the  magnetization  curve  is  quite  smooth  in 
this  region  (See  Fig.  4).  At  very  low  temperatures 
(1.4°K)  the  dhcp  form  also  shows  a  hysteresis  loop  which 


would  imoly  ferrimagnetic  or  ferromagnetic  ordering.  The 
effective  moment  calculated  from  the  susceptibility 

data  in  the  paramagnetic  region  for  the  two  forms  is  near 


68  which  compares  favorably  with  the  theoretical  value  of 
70.8  computed  using  g(j(J+l) )1//2  of  3-62  and  7*6  for  tri¬ 


positive  Nd  and  Tm  respectively. 


Lock1 7) 


investigated  the  magnetic  behavior  of  Sm- 


metal  and  noticed  a  sharp  decrease  in  the  susceptibility. 
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characteristic  of  antiferromagnetic  ordering  at  15°K,  and 
a  kink  in  the  curve  between  110°K  to  150°K.  Anomalies  in 
heat  capacity  have  been  reported  by  Jennings  et  al.1^  at 
105°  and  13.8°K.  Features  similar  to  those  observed  by 
Lock  were  obtained  by  Jayaraman  and  Sherwood  in  the  mag¬ 
netization  curve  for  the  normal  form  of  Sm.  For  the  pres¬ 
sure-induced  dhcp  form  there  was  a  much  more  pronounced 
kink  in  magnetization  around  100 °K.  Further ,  the  peak 
(See  Fig.  5)  occurring  at  the  lower  temperature  has  shifted 
to  27°K  and  is  somewhat  smaller  in  magnitude. 

iq^ 

McWhan  and  Stevens  have  suggested  that  the  double 
peaks  observed  in  the  A.C.  susceptibility  of  Sm-form  of 
Gd,  Tb  and  Dy  represent  ordering  of  hexagonal  (H)  layers 
and  cubic  (C)  layers  at  different  temperatures,  in  analogy 
with  the  behavior  of  Nd.  In  the  latter  metal  neutron  dif¬ 
fraction  studies20^  in  the  range  1.6°  to  20°K  have  shown 
that  the  hexagonal  site  layers  (H)  and  the  cubic  site 
layers  (C)  present  in  the  structure  order  at  different 
temperatures,  the  former  (H)  at  19°K  and  the  latter  (C)  at 
7.5°K.  McWhan,  Jayaraman  and  Sherwood  (to  be  published) 
suggest  that  in  the  Sm-metal  also,  the  cubic  and  hexagonal 
layers  order  at  different  temperatures.  The  kink  in  the 
magnetization  curve  at  the  higher  temperature  (100°K-3.20°K) 
we  believe  is  due  to  ordering  of  the  (C)  layers  and  the  low 
temperature  peak  (27°K  and  1.4°K)  associated  with  the  order¬ 
ing  of  (H)  layers.  The  specific  heat  measurements  of 
Jennings  et  al.  give  a  considerably  higher  entropy  under 
the  105°K  maximum  in  the  Cp  than  that  under  the  low  temp¬ 
erature  maximum,  which  would  seem  to  suggest  the  opposite 
sequence  for  the  ordering.  Our  suggestion  is  based  on  the 
following  facts.  In  the  normal  form  (rhombohedral  struc¬ 
ture)  the  hexagonal  to  cubic  layer  ratio  is  3  to  1  while 
in  the  dhcp  it  is  1  to  1.  Since  the  kink  in  the 
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magnetization  curve  at  the  higher  temperature  is  much  more 
pronounced  for  the  dhcp  form,  it  must  be  associated  with 
the  ordering  of  the  cubic  layers  for,  the  ratio  of  (C)  to 
(H)  layers  increases  in  going  from  the  normal  form  to  the 
dhcp  form  of  Sm.  It  would  also  follow  that  in  the  dhcp 
form  the  peak  associated  with  ordering  of  (H)  layers  should 
diminish.  This  is  precisely  what  is  observed. 

We  wish  to  thank  D.  B.  McWhan,  T.  H.  Geballe  and 
H.  J.  Williams  for  discussions  and  R.  G.  Maines  for  tech¬ 
nical  assistance  in  the  high  pressure  experiments. 
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VALENCE  BOND  FORMATION  IN  THE  CLOSE -PACKED 

RARE  EARTH  METALS 


Forrest  L.  Carter 

U.S.  Naval  Research  Laboratory,  Washington,  D>  C. 

ABSTRACT 

A  modified  Valence  Bond  treatment  of  the  close -packed  metals 
has  been  developed  subject  to  the  restrictions  that:  1)  all  near  neigh¬ 
bors  are  bound  simultaneously,  2)  equivalent  atoms  form  equivalent 
bonds,  and  3)  certain  periodic  boundary  conditions.  The  use  of  bidi¬ 
rectional  orbitals  having  the  general  characteristics  ap  +/T-a^cL 

X  /  XV 

✓ 

give  rise  to  various  'phase  structure..'  involving  'cycles'  of  old  VB 
theory.  In  the  ccp  and  hep  structures  this  leads  to  a  natural  division  of 
the  electrons  and  orbitals  into  those  associated  with  1)  a  band  of  co¬ 
valent  bonds,  2)  a  MO  band  having  s -character  for  free  electrons  and 
3)  localized  orbitals  not  directed  toward  near  neighbors  and  suitable 
for  containing  unshared  electrons. 

In  the  ccp  structure  the  pertinent  orthonormal  set  of  bonding 
orbitals  contain  50%  d  character  (a  =  y^-)  whereas  in  the  hep  the  p  and 
d  character  of  the  bidirectional  orbital  can  change,  giving  rise  to  a 
variation  of  C/A  with  d  character.  At  d  character  of  ca.  58%  the  C/A 
value  has  a  minimum  of  approximately  1.54  and  rises  to  the  ideal  C/A 
=  1.63  at  a  d  character  of  50%  and  67%. 

The  trivalent  rare  earth  metals  form  a  near  ideal  system  to 
test  these  ideas  inasmuch  as  the  local  states  of  s,  p,  or  d  character 
are  unoccupied  and  the  three  electrons  per  atom  appear  to  be  involved 
primarily  in  the  covalent  band  rather  than  a  free  electron  band. 

The  presence  of  both  hexagonal  and  cubic  closest  packing  in  the 
light  (Ce  — ?Sm)  rare  earth  metal  structures  is  in  accord  with  the  ob¬ 
served  C/A  of  ca.  1.62  and  percent  d  character  of  ca.  50%  in  the  bi¬ 
furcated  bonding  orbitals.  Thereafter  the  observed  C/A  drops  to  a 
minimum  near  Er  and  starts  to  rise  again  at  Lu.  The  general  behavior 
is  the  same  as  indicated  by  the  theory  for  increasing  amounts  of  d 
participation  in  the  bonding  orbitals.  The  latter  is  expected  to  occur 
with  increasing  atomic  number  because  the  d  energy  levels  fall  more 
rapidly  than  the  p  levels  and  are  accordingly  used  to  a  greater  extent 
in  (shared  electron)  orbital  occupation. 
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Introduction 

Valence  Bond  treatments  of  metals  occur  rather  infrequently 
and  are  becoming  less  popular,  with  the  possible  exception  of  Paul- 
ings's  semi -empirical  discussion  of  the  metallic  radii.  ^  The  reason 
for  the  decline  is  due  in  part  to  the  immense  computational  difficulties 
of  a  VB  treatment  of  a  solid.  However  the  decline  may  be  attributed 
primarily  to  the  success  of  band  treatments  which  not  only  are  compu¬ 
tationally  simpler  but  also  provide  an  energy  estimate  for  excited 
states;  information  of  considerable  import  for  those  interested  in 
transport  problems.  Nevertheless  the  VB  treatment  of  molecules  re¬ 
sults  in  ground  state  calculations  which  are  significantly  better  than 
the  MO  (or  'molecular'  band)  treatment  as  a  direct  consequence  of 
permitting  greater  correlation  between  electrons. 

Accordingly  we  might  renew  our  faith  in  the  directed  bond  of 
VB  theory  and  take  courage,  as  some  chemists^)  have, from  the  belief 
that  the  glue  that  holds  together  the  carbon  atoms  of  diamond  is  also 
responsible  for  cohesion  in  metals.  To  this  end  we  note  that  in  terms 
of  the  simpler  concepts  of  VB  theory  in  combination  with  the  bidirec¬ 
tional  bond  introduced 3)  previously  that  we  can  give  a  simple  quali¬ 
tative  explanation  of  the  observed  variation  of  C/A  in  the  hep  rare 
earth  metals.  In  short  the  variation  of  C/A  is  a  result  of  the  increas¬ 
ing  stability  of  the  5d  orbital  relative  to  the  6p  orbital  as  a  function  of 
increasing  atomic  number. 

Electron  Classes  -  Although  proper  quantum  theory  forbids  the  label¬ 
ing  or  identifying  of  individual  electrons  or  electron  groups,  it  is  a 
wide  spread  practice  and  conceptually  a  very  useful  one  to  talk  about 
the  "core  electrons"  or  the  "valence  electrons",  etc.  Accordingly  the 
total  energy  of  a  crystalline  solid  may  be  similarly  divided  as  in 
Eqn.  1. 


"Total 


E  E 

Rare  Gas  Core  +  Local  Unshared  Electrons 

(i) 

+  ^Bonds  +  ^Bands  +  ^Magnetic  +  ^"etc. 

(Free  Electrons) 


The  first  term  of  the  equation  is  obvious  and  needs  no  discus¬ 
sion.  The  second,  third,  and  fourth  terms,  however,  are  of  par¬ 
ticular  interest  in  that  they  are  intimately  interconnedted  (as  we  shall 
see)  via  the  orthogonality  conditions  occuring  in  the  formation  of  bi¬ 
directional  bonding  orbitals.  The  second  term  will,  in  the  rare  earth 
metals,  be  related  primarily  to  the  electrostatic  energy  of  unshared, 
paired  and/or  unpaired,  electrons  of  the  f  orbitals  because  the  two 
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localized  orbitals  of  d-p  character  will  not  be  occupied.  These  local¬ 
ized  orbitals  will  have  unfavorable  directionality  for  bonding  or  will 
have  small  radial  extension. 

The  bond  orbitals  of  the  present  discussion  are  bidirectional 
and  accordingly  have  lobes  directed  simultaneously  to  two  near  neigh- 
bors.  We  may  expect  crude  measures  of  strength,  such  as  melting 
point,  etc.  to  be  related  to  the  overlap  and  promotion  energy  of  such 
orbitals.  Another  feature  of  the  present  treatment  is  that  in  addition 
to  the  localized  orbitals  and  bonding  orbitals,  other  orbitals  largely 
of  s  character  remain.  With  these  we  are  free  to  construct  Bloch 
type  wave  functions  for  development  of  a  modified  free  electron 
treatment. 

The  E^agnetic)  term  is  considered  to  arise  from  the  inter¬ 
action  of  clouds  of  diffuse  magnetic  dipole  distribution  of  primarily 
localized  electrons  and  not  to  the  electron  repulsion  effects  that  are 
the  source  of  the  magnetic  dipole.  While  these  terms  are  of  consider¬ 
able  interest  they  will  not  be  discussed  further  here. 

It  is  clear  that  the  fourteen  rare  earth  metals  are  going  to  con¬ 
stitute  a  very  important  testing  ground  for  numerous  physical  theories. 
Further  we  believe  that  this  series  constitutes  a  unique  test  for  a 
theory  of  chemical  bonds  in  metals.  That  uniqueness  resides  in  the 
following  features  of  the  rare  earth  metals:  1.  In  general  the  three 
valence  electrons  are  only  weakly  interacting  with  the  unshared  elec¬ 
trons,  a  situation  which  is  not  true  in  the  transition  metals.  2.  A 
gradual  change  takes  place  in  the  relative  energies  of  atomic  levels 
involving  5d,  6s,  and  6p  states.  3.  Most  important,  the  three  valence 
electrons  (excluding  Ce,  Eu)  are  primarily  involved  in  bond  formation, 
in  contradistinction  to  being  involved  in  a  free  electron  type  role.  In 
support  of  the  last  item  we  note  the  following:  (a)  The  rare  earth 
metals  are  metallic  conductors  of  unusually  high  resistivity.4)  Speak¬ 
ing  loosely  we  might  say  that  the  conducting  electrons  are  very  much 
involved  in  the  local  picture  (or  are  frequently  scattered  compared  to 
the  free  electron  treatment),  (b)  Augmented  plane  wave  calculations 
are  in  distinct  disagreement  with  the  free  electron  model.  ^)  (c)  Heats 
of  atomization  for  the  elements  at  STP  are  fully  consistent  with  the  in¬ 
volvement  of  three  electrons  per  atom  in  strong  bond  formation.2) 

Bonding  in  the  Close -Packed  Metal  Structures  -  Bidirectional  bonds 
of  the  same  form  can  be  used  to  discuss  the  bonding  in  both  the  cubic 
and  hexagonal  structures.  This  discussion  presumes  that  the  6s,  6p, 
and  5d  orbitals  are  sufficiently  close  in  energy  to  permit  bond  hybrid¬ 
ization,  however,  let  us  defer  the  consideration  of  f  orbitals  till 
later. 
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Fig.  1.  A  C-orbital  of  the  form  l//2  p  +  l//2  d  is  shown  in  the  (001) 
plane  of  the  fee  structure.  In  part  b  we  see  howx^uch  C-orbitals  on  near¬ 
by  atoms  can  be  used  to  form  a  cycle  involving  four  atoms  with  good 
overlap  between  them. 

Cubic  close -packed  bonding.  In  the  ccp  structures  each  atom 
has  12  near  neighbors,  with  four  of  these  forming  a  square  array  in 
each  principle  plane.  This  is  illustrated  in  Fig.  la  where  a  bidirec¬ 
tional  bond  of  the  form  apx  +  )f  1-a^  is  used  to  bond  the  two  atoms 

to  the  right.  If  a  =  a  second  orbital  equivalent  and  orthogonal  to 

the  first  may  be  used  to  form  bonds  to  the  two  atoms  on  the  left.  Sim¬ 
ilar  and  mutually  orthogonal  bonds  may  be  formed  in  the  two  other 
perpendicular  planes  for  bonding  the  other  eight  near  neighbors.  Such 
orbitals  correspond  to  an  electron  density  distribution  of  C^v  sym¬ 
metry  and  have  been  indicated  by  the  term  C-orbital  in  a  discussion  of 
the  rare  earth  chalcogenides  having  the  TI13P4  structure. The  six 
C-orbitals  of  the  ccp  structure  have  been  also  used  in  a  consideration 
of  bonding  in  the  distorted  icosahedral  coordination  of  the  (3 -tungsten 
structure. 


If  the  six  C-orbitals  comprise  the  bonding  orbitals  employing 
the  d  ,  d  ,  d  ,  px,  p  ,  and  p  atomic  orbitals  in  hybrid  formation, 
then’tXe  remaining  orbitXls  are  tfie  d^??  t  &z2>  and  s  states. 


The  d^y2 


and  d  ? 


orbitals  are  local  in  character  in  that  they 


are  directed  away  from  all  near  neighbors.  The  s  orbital  is  of  course 
non -directional  and  may  be  used  as  a  basis  for  a  Bloch  function  which, 
at  the  atomic  sites,  is  orthogonal  to  the  bonding  orbitals  in  the  first 
approximation. 


Hexagonal  close -packed  bonding.  The  same  general  situation 
obtains  in  the  hexagonal  close  packed  structures.  Six  C-type  orbitals 
(somewhat  distorted  from  C2v  symmetry)  are  used  to  bond  the  twelve 
near  neighbors,  two  per  C-orbital,  one  neighbor  in  the  basal  plane  and 
one  neighbor  above  or  below  the  basal  plane.  These  orbitals  cannot 
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Fig.  2.  This  shows  the  12  near¬ 
est  neighbors  of  the  central  atom 
in  the  hep  arrangement;  6  in  the 
basal  plane,  three  above,  and 
three  below.  The  C-orbitals 
are  represented  schematically 
by  the  shape  of  their  large 
lobes  in  the  plane  containing 
the  central  atom  and  the  lobe 
maxima;  the  small  lobes  are 
not  indicated.  Each  of  the  six 
orthogonal  C-orbitals  binds 
two  atoms  which  are  at  right 
angles  at  the  central  atom, 
e.  g.  No.  1  and  No.  1'. 


be  obtained  by  inspection  however,  as  can  those  of  the  cubic  structure, 
and  we  must  use  more  general  methods.  Let  us  take  a  function  of  the 
form  Eqn.  2  and  from  it 

©  =  A  s  +  Bp  +  Cp  +  Dp  +  Ed  +  Fd  +  Gd 

x  y  z  xy  xz  yz 

+  Hd^_y2  +  Idz2  (2) 

obtain  five  others  according  to  the  point  symmetry  32.  If  these  six 
general  orbitals  are  to  be  mutually  orthogonal  and  normalized  the 
above  coefficients  must  satisfy  the  following  conditions: 


A2  +  I2  =  D2  -  1/6  (3) 

B2  +  (j  +  H2  =  1/3  (4) 

C2  +  E2  +  F2  =  1/3  (5) 

BC  =  EH  +  FG  ‘  (6) 


Inasmuch  as  the  3  fold  axis  is  taken  along  the  Z  direction  we 
see  that  except  for  sign  the  pz  component  will  be  the  same  for  each  of 
the  6  orbitals  (the  same  will  be  true  for  the  s  and  dz2  orbitals).  Ac¬ 
cordingly  the  pz  orbital  will  be  employed  in  its  entirety  in  the  forma¬ 
tion  of  the  orbitals  corresponding  to  the  point  group  32.  At  our  cur¬ 
rent  level  of  approximation  the  same  will  be  roughly  true  of  the  dz2 
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Fig.  3.  This  indicates 
the  shape  of  a  C-orbital 
in  the  plane  of  its  lobe 
maxima.  This  orbital 
contains  66.  6%  d  char¬ 
acter  and  can  be  used 
for  bonding  in  the  hep 
structure.  The  coef¬ 
ficients,  in  terms  of 
Eqn.  2,  are  A=0,  B=-C 
=-F=-G=l//12,  JD=E=-H-I 
•-1//6. 


orbital.  This  follows  from  Eqn.  3  and  the  fact  that  s  character  in  the 
bond  will  be  small  because  it  is  spherical  and  has  small  directed  ra¬ 
dial  extension;  therefore  overlap  will  fall  off  with  large  amounts  of  s 
character. 

In  Fig.  2  we  have  schematically  illustrated  the  relationship  of 
the  various  C -orbitals  of  Eqns.  2  to  6  to  one  another  and  the  twelve 
nearest  neighbors  in  the  hep  structure.  For  example,  atom  No.  1  of 
the  basal  plane  and  atom  No.  1'  of  the  upper  three  atoms  are  joined  by 
a  C-orbital.  This  orbital  is  related  to  the  orbital  joining  atoms  No.  2 
and  2'  by  a  120°  rotation  about  the  3  fold  axis  (z  direction).  The  or¬ 
bitals  joining  atoms  No.  2  and  2'  and  No.  3  and  3’  are  similarly  re¬ 
lated.  The  C-orbital  whose  lobes  are  directed  toward  atoms  No.  4 
and  4'  is  related  to  C-orbital  1-1'  by  a  2  fold  rotation  about  the  x  axis. 
The  two  remaining  orbitals  not  shown  are  related  to  C-orbital  4-4*  by 
the  3  fold  rotation  axis. 

If  we  inspect  a  C-orbital  of  the  hep  structure  in  the  plane  con¬ 
taining  the  origin  and  the  two  maxima  of  the  lobes  as  in  Fig.  3  we  see 
that  indeed  it  is  very  similar  to  the  C-orbital  of  the  ccp  structures, 
Fig.  1.  There  is  another  important  point  to  note  in  Fig.  3,  i.  e.  ,  the 
atoms  joined  by  the  C-orbital  form  a  right  angle;  hence  this  is  true  in 
loth  the  hep  and  ccp  structures. 

Having  formed  six  C -orbitals  for  bonding,  from  a  basis  of  nine 
atomic  orbitals,  it  is  clear  that  three  orbitals  remain.  One  of  these, 
composed  of  6s  plus  5dz2  character,  will  be  suitable  for  the  atomic 
part  of  a  Bloch  function.  The  d  2  character  here  will  be  small,  of 
course,  if  little  s  character  is  used  in  the  C -orbitals.  The  other  two 
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orbitals,  of  5d  and  6p  character,  may  be  expected  to  have  small  over¬ 
lap  with  neighboring  atoms  and  hence  can  be  described  as  local  in 
character. 

We  have  argued  above  on  the  basis  of  electrical  conductivity, 
heats  of  atomization,  etc.  ,  that  small  occupation  is  expected  of  a 
band  (composed  here  of  6s  +  5dz2  character)  having  a  free  electron 
aspect.  In  addition  there  is  even  less  reason  to  suppose  that  the-.lo- 
calized  orbitals  would  be  occupied. 

We  may  then  suggest  that  since  the  bonding  electrons  of  the 
heavier  rare  earths  will  have  small  interaction  with  either  the  loca¬ 
lized  electrons  or  the  Bloch  electrons  these  metals  constitute  an  ex¬ 
cellent  testing  ground  for  the  theories  of  bonding  in  metals. 

Variation  of  Bond  Strength  with  d  character  -  Earlier  it  was  noted 
that  the  two  atoms  bound  by  the  C -orbitals  form  a  right  angle.  Using 
a  generalized  notion**)  of  Pauling's  concept  of  the  bond  'strength'  S  of 
an  orbital  we  may  examine  the  effect  of  varying  the  amount  of  d 
character  of  a  C-orbital.  This  amounts  to  taking  the  maximum  pro¬ 
jection  of  the  angular  component  of  the  lobe  on  the  pertinent  inter- 
nuclear  axis  as  is  indicated  in  Fig.  3.  For  an  oi*bital  of  pure 
symmetry  this  variation  is  illustrated  in  Fig.  4.  We  note  that  only 
for  small  and  large  amounts  of  d  character  does  the  strength  S  fall 
off  drastically.  Therefore  it  remains  to  determine  the  direction  of 
the  non -basal  plane  lobe  with  variation  in  d  character. 


no 


Variations  of  C/A  with  d-character  -  An  approximation  to  the  vari 
ation  of  the  hexagonal  C/A  ratio  with  the  amount  of  d  character  in  the 
C -orbital  is  obtainable  from  the  following  stratagem.  An  orbital  of 
the  type  ap1  +  yi-a^  d'  is  rotated  about  the  z1  axis  by  an  angle  oo 

(to  =  -45  )  until  the  maximum  of  the  positive  lobe  is  a  little  past  the  x* 
axis.  The  negative  large  lobe  then  lies  along  the  +y*  direction  (see 
Fig.  3).  The  orbital  is  then  rotated  about  the  x1  axis  by  an  angle  2T 
which  then  directs  the  negative  lobe  toward  one  of  the  atoms  above  the 


basal  plane. 


tan  2T~ 


The  angle  yis  related  to  the  C/A  ratio  by  Eqn.  7. 

-  (°  \  an 

"  va~7  2 


This  C -orbital  is  then  in  the  same  relationship  to  the  near  neighbors 
as  is  C-orbital  1-1'  of  Fig.  2  if  we  take  the  x'  axis  as  the  internuclear 
axis  between  the  central  atom  and  atom  No.  1  and  the  Z  and  Z1  axes 


as  the  same.  If  the  coefficients  of  p1  and  d'  ^  are  expressed  in  terms 
of  the  angles  oj  and  }T  and  set  equal  toz  z  then  to  can  be  elimi¬ 

nated  to  give  an  expression,  Eqn.  8,  between  y  and  a^,  the  amount  of 
p  character.  Z  f 

sin2  r  =  4“  \r~2  +  “X  )  (8) 


By  combining  Eqns.  7  and  8  we  plot  C/A  as  a  function  of  d 
character  in  Fig.  5.  Three  points  are  important:  (1)  a  given  C/A  does 
not  correspond  to  a  unique  value  of  d  character  in  the  C-orbital;  (2) 
large  amounts  of  p  character  (  ySQ%)  and  small  amounts  of  p  character 
(<  33%)  correspond  to  high  C/A  ratios  and  are  not  pertinent  to  a  C- 
orbital  description  since  the  lobes  become  quite  disimilar;  (3)  a  mini¬ 
mum  C/A  value  is  observed  at  58%  d  character.  C-orbitals  satisfying 
all  the  orthogonality  conditions  (Eqns.  3-6)  do  not  rise  as  fast  as  does 
the  curve  of  Fig.  5  near  high  values  of  d  character.  Similar  curves 
will  be  presented  at  the  conference  for  various  ratios  of  the  projected 
maxima  of  the  two  lobes  according  to  bond  order -bond  length  con- 
side  rations. 

C/A  in  the  Rare  Earth  Metals  -  As  the  atomic  number  of  the  rare 
earth  metal  increases  we  may  expect  the  5d  orbitals  to  become  increas 
ingly  stable  relative  to  the  6s  orbital  while  the  6p  state  will  be  relative 
ly  constant;  this  is  schematically  illustrated  in  Fig.  6.  In  view  of  our 
arguments  concerning  the  (non)  use  of  the  6s  orbital  in  bond  formation, 
we  are  really  interested  in  the  increasing  separation  between  6p  and 
5d  levels.  As  this  separation  increases  with  atomic  number  more  d 
character  will  be  involved  in  the  C-orbital  bonding  at  the  expense  of  p 
character.  According  to  Fig.  5  we  may  expect  the  C/A  ratio  to 
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decrease  to  a  minimum  and  then  to  increase  again  as  the  d  character 
approaches  67%.  Such  a  variation  of  C/A  is  observed  in  the  heavier 
trivalent  rare  earths  (see  Fig.  7). 

Another  point  of  interest  is  that  at  50%  d  character  Fig.  5 
shows  the  C/A  to  be  ideal  (i.  e.  ,  1.63)  which  is  the  same  value  as  occurs 
in  the  ccp  structure.  However  since  orthogonality  conditions  require 
the  C-orbital  to  also  have  50%  d  character  in  the  cubic  structures  we 
have  a  continuous  transition  from  the  ccp  to  hep  structures  in  terms 
of  the  composition  of  the  C-orbital.  This  is  essentially  observed  in 
the  layer  phases  of  the  lighter  rare  earths  (Fig.  7)  which  contain  both 
cubic  and  hexagonal  close  packing  with  a  C/A  just  a  little  below  the 
ideal  value  of  1.63.  In  the  heavier  rare  earth  metals  we  may  expect 
less  participation  in  layer  structures  because  then  a  larger  promotion 
energy  would  be  associated  with  the  cubic  sites. 

F  Orbital  Participation  in  Bonding  -  In  Fig.  6  we  note  that  the  ener¬ 
gies  of  the  6p,  6s,  5d,  and  4f  orbitals  are  rather  close  in  the  lighter 
rare  earths.  It  would  not  be  surprising  then  if  f  orbitals  played  a  role 
in  the  bonding.  In  particular  we  remark  here  that  the  f  orbitals  are 
of  odd  symmetry  like  the  p  orbitals  and  can  be  employed  in  a  similar 
fashion  in  the  formation  of  C-orbitals.  Infact  they  can  be  used  in  com¬ 
bination  with  cL.  type  orbitals  to  give  lobes  whose  maxima  make  an 
angle  of  90  or  somewhat  less,  whereas  the  p  orbital  does  not  have 
this  property.  In  this  possibility  may  reside  the  explanation  for  the 
C/Anot  declining  with  atomic  number  before  Sm.  In  the  heavier  rare 
earths  however  we  may  presume  that  the  radial  part  of  the  f  orbital  is 
so  concentrated  that  it  no  longer  is  able  to  make  a  contribution  to  the 
bonding. 

Final  Remarks  -  The  bidirectional  nature  of  the  C-orbital  permits  the 
atom  in  question  to  bind  two  neighbors  per  C-orbital  if  the  signs  (or 
phases)  of  the  overlapping  lobes  are  the  same.  By  requiring  the  match¬ 
ing  of  phases,  atoms  are  linked  together  in  chains  or  rings  of  even 
number  of  atoms.  Rings  or  cycles  involving  four  atoms  could  be  form¬ 
ed  in  the  ccp  structures  (Fig.  lb)  and  would  lead  to  strong  bond  form  ¬ 
ation.  Such  rings  can  also  be  formed  in  the  hep  structures  but  appear 
to  have  somewhat  less  overlap.  Cycles,  as  rings  or  chains,  can  be 
used  to  form 'phase  structures' ',  in  which  all  atoms  are  consis¬ 
tently  bonded  and  which  may  or  may  not  have  the  same  periodicity  of 
the  lattice. 

It  is  necessary  to  say  something  about  the  symmetry  of  the  C- 
orbitals  in  relation  to  the  higher  symmetry  we  expect  for  electron 
bond  density.  We  note  that  the  C-orbitals  in  Fig.  2  are  related  by 
point  symmetry  32  whereas  the  atom  site  in  the  crystal  has  higher 
symmetry  due  to  the  presence  of  mirror  planes.  This  symmetry  is 
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restored  if  we  consider  the  mirror  image  of  Fig.  2  is  equally  likely. 
Phase  structures  containing  one  or  both  enantimorphic  forms  of  Fig.  2 
can  be  devised.  The  same  situation  applies  to  the  C -orbitals  in  ccp 
positions.  It  is  worth  noting  further  that  transitions  between  hep  and 
ccp  layers  can  be  discussed  by  means  of  cycles  and  phase  structures 
where  each  atom  is  consistently  bonded  (in  terms  of  phase  signs  and 
orientation  of  lobes)  to  all  its  neighbors. 

Finally  we  note  that  the  use  of  this  form  of  covalent  binding 
does  not  deny  these  metals  their  electrical  conductivity  since  the 
number  of  electrons  per  atom  is  insufficient  to  fill  these  'covalent 
bands'  ^). 
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